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ABSTRACT : 



An objective-numerical scheme for locating and specifying the inten- 
sity of hyperbaroc linic zones (HBZs) is presented and applied in both 
research and operational environments. The operational model, called 
DG0*, employs the isobaric gradient of virtual potential temperature (0*) 
to specify baroc linic ity while a directional second derivative of 9* is 
the parameter used to locate the HBZs and indirectly give a measure of 
their intensity. 

The DG0* model is diagnosed from an analysis of modeled temperature 
profiles and a number of real-atmosphere temperature fields at 1000, 850, 
700, and 500 mb, as taken from winter and summer situations in the years 
1964-67. Consideration is given to location, intensity, frequency, slope, 
and vertical and horizontal distribution and consistency of the HBZs in 
relation to the numerical processing of data. A key synoptic time, 0000 
GMT 6 March 1966, serves to exemplify analytical and diagnostic features 
of the HBZs. The numerically-analyzed 1000 mb fronts, produced operationally 
by Fleet Numerical Weather Central, Monterey, California since 1965, are 
statistically and graphically compared to subjectively-analyzed versions 
for winter and summer periods up to 30 days. The relation of the HBZ to 
vertical motion and cloud is illustrated by a case study from December 1966. 

Experiments, using a high-resolution scheme for analyzing the HBZs, are 
also presented. A proposed model, employing parameters of vertical sta- 
bility and isentropic slope, utilizes the thermal detail inherent in 
vertically-resolved radiosonde data to give an enhanced field of baro- 
clinicity on the constant pressure surface. In addition, it is shown how 
the model may be coupled with a vertical extrapolation procedure to at 
least triple the amount of baroclinicity data used for HBZ analysis. 



ACKNOWLEDGMENTS 



Space hardly permits citing all of the individuals and groups 
who helped to bring the research to this terminal point. Without the 
technical and financial support of the Navy’s Fleet Numerical Weather 
Central, Monterey, California (FNWC) , and its leader. Captain Paul 
Wolff, the study would not have begun. Mr. Leo Clarke, research 
meteorologist at FNWC, wrote most of the computer programs as well as 
collaborated in ensuring utility of the model for the objective analysis 
of fronts by the practicing meteorologist. In addition, much impetus 
in the early stages was given by financial support from the Office of 
Naval Research; more recently the project was supported by the Navy 
Weather Research Facility, Norfolk, Virginia. Appreciation is expressed 
to the Naval Postgraduate School (PGS) for the facilities afforded this 
investigator as a member of the faculty. In particular. Dr. George J. 
Haltiner, Chairman of the Department of Meteorology, PGS, is thanked for 
his constant encouragement to pursue the goals of this investigation. 

For many of those facets of the problem that involved feasibility 
experimentation and model testing in both time and space scales, many, 
many students in the Department of Meteorology, PGS, gave valuable 
assistance, as did the civilian and military technicians associated 
with the teaching and research programs at the School. 

Last, but certainly not least. Dr. Noel E. LaSeur , as well as 
other members of the Department of Meteorology, Florida State University, 
are sincerely thanked for giving of their time and advice in fostering a 
research study which was accepted by Florida State as a suitable doctoral 
dissertation study. 

ii 



TABLE OF CONTENTS 



Page 



ABSTRACT i 

ACKNOWLEDGMENTS ii 

LIST OF TABLES vi 

LIST OF ILLUSTRATIONS viii 

Chapter 

I. HISTORICAL BACKGROUND AND OBJECTIVES OF STUDY 1 

II. HYPERBAROCLINIC-ZONE ANALYSIS: DEVELOPMENT AND 

CHARACTERISTICS OF A NUMERICAL MODEL 9 

The Mode 1 



A One “Dimensional Schematic Analysis of DG9* 

Analysis Characteristics of the Model 
The Simulation Experiments 
Frontal-Zone Intensity 
Variation in Mesh Length 
Closely Adjacent Fronts 
Vertical Slope of HBZ Boundaries 
Use of Layer-Mean Temperature, 0*, in the DG0* 

Ana lysis 

General Conclusions 

III. AN OPERATIONAL APPLICATION OF THE HYPERBAROCLINIC-ZONE 

ANALYSIS MODEL 36 

Data Processing Framework for the Operational 
DG0* Analysis 

An Example of Two-Dimensional Analyses of DG0* 

Operational Hemispheric HBZ Analysis 

Objective Frontal Analysis for 0000 GMT 6 March 1966 
Objective Frontal Prognosis for 1200 GMT 7 March 1966 
Comparison of Subjective and Numerical Frontal Analyses 
1000 mb Analyses 
850 mb Analyses 

Modification of DG0* to Enhance the An^alyses of Occlusions 
Effect of Mesh-Length Variation on DG0* 

Summary 



iii 



IV. 



EXPERIMENTS TO ENHANCE THE RESOLUTION OF 
BAROCLINICITY 



Page 



71 



Computation of G0 and DG0 

The Computation of SO/Sp 
The Computation of 6p/6n 
G0* and its Comparison to G0* 

DG0* and its Comparison to DG0* 

Recomputation of G0 with Certain Approximations 
Removed 

Cross Section Analyses 
G0 Extrapolations 

Analysis of Combined G0* and G0** Data 

c c 

Summary 

V. THE DG0* MODEL: CASE STUDIES OF THE SPACE AND TIME 

CONTINUITY OF HYPERBAROCLINIC ZONES 107 

A Summer-Season Case Study 

DG0^" Analyses at 1000, 850, 700, and 500 mb 
Vertical Continuity of HBZ Boundaries 
A Winter-Season Case Study 

HBZ Boundaries in Vertical Cross Section 

VI. THE DG0* MODEL: A CASE STUDY OF HYPERBAROCLINIC ZONES 

IN REIATION TO CLOUD AND VERTICAL MOTION 132 

Selection of Period and Area of Study 
Patterns of Vertical Motion and Cloud in Relation 
to the Numerically-Analyzed Fronts 
Satellite-Observed Cloud Bands Versus Hy perbaroc linic 
Zones 

Stratification by Latitude 
Stratification by Geographical Area 

VII. THE DG0* MODEL: HEMISPHERIC FREQUENCY AND AVERAGE 

INTENSITY OF HYPERBAROCLINIC ZONES ..... 149 

Winter-Season Sample: 15 January through 

14 February 1966 

Numerical Frequency 
Numerical Intensity 
NMC Frequency 

Summer-Season Sample: July 1967 

Numerical Frequency 
Numerical Intensity 
NMC Frequency 



iv 



Page 



VIII. RESUME’ AND PROPOSALS FOR FURTHER STUDY 163 

Summary of Results 

Proposed Areas for Further Research 
Final Remarks 

APPENDIX A - PRELIMINARY EXPERIMENTS IN HYPERBAROC LING -ZONE 

ANALYSIS coco 169 

Experiments Using Wet-Bulb Potential Temperature 
Relations of 0 and V0 to the Polar Front 
Summer Case §tudy: Manual Analysis 

Winter Case Study: Numerical Analysis 

Experiments Using Equivalent Potential Temperature 
Summary 

REFERENCES ‘ 184 

INITIAL DISTRIBUTION LIST 190 

FORM DD 1473 I 95 



V 



LIST OF TABLES 



Tab le Page 

1-5 Grid Point Values of 0* and Computed G0* and 

DG0* for Modeled 9* Profiles of Fig. 3 

(1) N, M, and W Profiles; HBZ A0* = IOC; 3- and 

5-Point Differentiations 16 

(2) N Profile; HBZ A0* = 5C ; 3- and 5-Point 

Differentiations 20 

(3) N and M Profiles; Mesh Length Halved; HBZ 

A0* = IOC; 5-Point Differentiation 23 

(4) Horizontally Combined N and M Profiles for 
Multiple Frontal-Zone Cases; HBZ A0* = IOC; 

5-Point Differentiation 25 

(5) Vertically Combined N, M, and W Profiles for 

Mean 1000/700 mb 0* and Frontal-Zone Slopes 
of 1/120 and 1/240; HBZ A0* = IOC; 5-Point 
Differentiation 30 

6 Statistical Relations Between Numerically- 
Analyzed 1000 mb and Manually-Analyzed NMC 
Surface Fronts for Period 23 August--6 September 

1965 57 

7 850 mb AG0, 0000 GMT 6 March 1966 78 

8 850 mb AG0’ , 0000 GMT 6 March 1966 92 

9 850 mb 6p/6n and G0 Computations, 0000 GMT 

6 March 1966 98 

10 Space-Continuity Statistics for Numerically- 

Analyzed HBZs in Cross Sections of Figs. 32-34 . . . 126 

11 Time-Continuity Statistics for Numerically- 

Analyzed HBZs in Cross Sections of Figs. 32-34 . . . 129 

12 Relation of Satellite-Observed Clouds to 
Hyperbaroclinic and Subbaroc linic Zones, 

7-17 December 1966 145 

vi 



Table 



Page 



13 


Relation of Satellite-Observed Cloud Bands to 
Hyperbaroclinic and Subbar oc linic Zones, 7-17 
December 1966 


147 


14 


Monthly Mean 0 and 0 and their Standard 
we 

Deviations for the Polar Front 


170 


15 


Percentage Frequency of 0^ Along the 850 mb 

Polar Fronts Analyzed by NMC and CAO 


172 


16 


Average 0^, as a Function of Longitude, Along 

the 850 mb Polar Fronts Analyzed by NMC and CAO . . 


173 


17 


Table of Related Values from Eqs , (27) and (29), 
as Applied to a February 850 mb Polar Front .... 


179 



VLl 



LIST OF ILLUSTRATIONS 



Figure Page 

1 Surface Fronts for 0000 GMT 5 March 1964, as 
Analyzed by Sixteen National and International 

Weather Centers 6 

2 Schematic Profiles of 9*, G0*, and DG0* 11 

3 Modeled Profiles of 0* 14 

4 Modeled and Numerically-Computed HBZ Boundaries 

in Vertical Cross Section 28 

5 Typical Plan Views and Profiles of 0*, G0*, and 

DG0* 40 

6 1000/700 mb DG0* Analysis, 0000 GMT 6 March 1966, 

with Superimposed NMC Surface Fronts 43 

7 NMC 1000/500 mb Thickness Pattern, 0000 GMT 6 March 

1966, with Superimposed NMC Surface Fronts .... 46 

8 Operational 1000/700 mb DG0* Analysis, 0000 GMT 

6 March 1966 49 

9 FNWC Sea-Level Pressure Field, 0000 GMT 6 March 

1966, and the Numerically-Analyzed 1000 mb Fronts , 51 

10 FNWC 850 mb Isotherm and Contour Fields, 0000 GMT 

6 March 1966, and the Numerically-Analyzed 1000 mb 
Fronts 53 

11 36-Hour 1000 mb DG0* Prognosis and Verifying 

Analysis, 1200 GMT 7 March 1966 55 

12 850 mb DG0* Analysis, 0000 GMT 6 March 1966 .... 59 

13 1000/700 mb GK0* Analysis, 0000 GMT 6 March 1966, 

with Superimposed NMC Surface Fronts 61 

14 1000/700 mb (DG0* + GK0*) Analysis, 0000 GMT 6 

March 1966, with Superimposed NMC Surface Fronts . . 63 



viii 



Figure Page 

15 Modification to 1000/700 mb DG0^ Analysis for 

Frontal Occlusion, 0000 GMT 6 March 1966 64 

16 1000/700 mb DG0* Analysis over Selected Areas 
as a Function of Mesh Length, 0000 GMT 6 March 

1966 67 

17 850 mb G0* and Subjectively-Analyzed G0’ , 0000 

GMT 6 March 1966 76 

18 850 mb Objectively-Analyzed G0*, 0000 GMT 6 March 

1966 82 

19 850 mb DG0* and DG0' , 0000 GMT 6 March 1966 83 

20 Various Subjective and Objective 850 mb Frontal 

Analyses, 0000 GMT 6 March 1966 85 

21 850 and 700 mb HBZ Boundaries from DG0* and DG0’ 

Analyses, 0000 GMT 6 March 1966 86 

22 850 mb G0* , 0000 GMT 6 March 1966 91 

c 

23 Vertical Cross Section of 0* for Stations 78367 

to 72655, 0000 GMT 6 March 1966 94 

24 Vertical Cross Section of 0 for Stations 78367 

to 72655, 0000 GMT 6 March 1966 95 

25 850 mb G0* and G0”, 0000 GMT 6 March 1966 102 

26 1000 mb DG0*, 0000 GMT 20 August 1965 110 

27 850 mb DG0*, 0000 GMT 20 August 1965 Ill 

28 700 mb DG0*, 0000 GMT 20 August 1965 113 

29 500 mb DG0*, 0000 GMT 20 August 1965 114 

30 1000, 850, 700, 500, and 300 mb Positive DG0^ 

Axes, 0000 GMT 20 August 1965 115 

31 Schematic 1000 mb Frontal Contour Showing Locations 

of Cross Section Groups in Figs. 32-34 118 

32 Vertical Cross Sections of HBZ Boundaries, Groups 

I and II 120 

Vertical Cross Sections of HBZ Boundaries, Group III . 122 

ix 



33 



Figure Page 

34 Vertical Cross Sections of HBZ Boundaries, 

Groups IV and V 124 

35 Schematic Vertical Cross Section of Coexisting 

Low- and High-Tropos phere HBZs 128 

36 Time Continuity of 850 mb Numerically-Analyzed 
Fronts 5 0000 GMT 26 January to 1200 GMT 27 January 

1967 131 

37 Schematic Illustration of Idealized Satellite- 
Observed Cloud Formations in Relation to an 

Occluded Surface Front 134 

38 Satellite-Observed Cloud Patterns, Vertical 
Motion Fields and Positive DG0* Axes , North 

Pacific Region, 15 December 1966 138 

39 Satellite-Observed Cloud Patterns, Vertical 
Motion Fields and Positive DG0* Axes , North 

Pacific Region, 16 December 1966 140 

40 Satellite-Observed Cloud Patterns, Vertical 
Motion Fields and Positive DG0* Axes, North 

Pacific Region, 17 December 1966 142 

41 Frequency and Intensity of Numerica lly -Analyzed 
1000 mb Fronts for the Period 15 January 

through 14 February 1966 152 

42 Frequency of Manually -Analyzed NMC Surface 
Fronts for the Period 15 January through 

14 February 1966 156 

43 Frequency and Intensity of Numerically- 

Analyzed 1000 mb Fronts for July 1967 158 

44 Frequency of Manually^Ana lyzed NMC Surface 

Fronts for July 1967 161 

45 850 mb Polar Front Positions, 1200 GMT 21 June 

1963: Various Subjective and Statistical Analyses . 175 

46 Portion of 850 mb 0* and 0* Analyses over the 

United States, 1200^GMT 5 February 1964 180 

47 Frontal Field, F*, Generated for Use in Locating 

the 850 mb Arctic Front, 0000 GMT 5 March 1964 . . . 182 



X 



NPS-51RD0031A 
23 March 1970 

HYPERBAROCLINIC ZONES: NUMERICAL ANALYSIS 

AND APPLICATION TO FRONTS* 

R. J. Renard 
CHAPTER I 

HISTORICAL BACKGROUND AND OBJECTIVES OF STUDY 

It is difficult to ascertain exactly when meteorologists first 
became aware of, or concerned with, atmospheric baroclinicity . According 
to Bergeron (1959), Kbppen was one of the first nineteenth century 
meteorologists to perceive observational evidence of air masses and their 
baroclinic boundaries, conventionally called fronts . However, it was not 
until the 1920s and 1930s that the so-called Norwegian air-mass and frontal 
concepts (Bjerknes, 1918) were firmly established. These ideas live on 
today in the form of frontal analyses (e.g., warm, cold and occluded 
fronts) found on operational weather charts of every major meteorological 
service in both hemispheres, without exception. 

Since the 1940s, the broader concept of a baroclinic frontal zone , 
vice a frontal line , has become an increasingly important theoretical and 
synoptic consideration in the life history of many important facets of 
weather systems. In this connection, the term hyperbaroc linic zone (HBZ) , 
after Godson (1951), is most appropriate to describe a restricted area of 
relatively large values of baroclinicity, without necessarily connoting a 
zone separating identifiable air masses. Many investigations in the pre- 
computer era have demonstrated the importance of HBZs to the dynamic and 
thermodynamic processes and structure of the atmosphere (Bjerknes, 1951; 
McClain and Danielsen, 1955; Newton, 1958; Palmen, 1951; Palmen and Newton, 
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1951; Reed and Sanders, 1953; Sanders, 1955). A synthesis of theories 
and examples of development and motion of waves, vortices, jet streams, and 
cloud and precipitation regimes as a function of baroc linicity may be found 
in dynamics texts (e.g., Haltiner and Martin [1957] or Hess [1959]) and 
applied manuals (e.g., Riehl et al. [1952, 1954], the United States Navy 
[1952] and the United States Weather Bureau [I960]). However, quanti- 
tative knowledge of the behavioral characteristics of HBZs on a synoptic 
hemispheric scale has remained incomplete. This is due to the short- 
comings imposed by subjective analysis with its attendant ambiguity in 
locating and depicting the important HBZs and to the derived, not observed, 
nature of baroc linicity . The operational convention of specifying only 
the warm-side boundary of the HBZ (Taljaard et al., 1961) has further 
detracted from attention to the zone rather than the line . 

Attempts at overcoming the frustrations of the manual analysis of 
HBZ boundaries have been pioneered by the Meteorological Service of Canada 
(Anderson et al., 1955). Their frontal-contour model, comprising three 
fronts and four air masses has been used to assist in a quasi-objective 
approach to frontal analysis. However, the model places undue emphasis 
on the labeling of air masses and fronts by means of empirically-determined 
frontal values of various temperature parameters, notably wet-bulb potential 
temperature. Moreover, the time involved in the subjective manipulations 
necessary for the objective placement of frontal contours exceeds the human 
capacity for anything but a regional effort, if operational deadlines are 
to be met. 

The dilemna of real-time urgency in frontal analysis suggests 
numerical-objective vice manua 1-sub jective or manua 1-objective analysis 
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of HBZs . Such an approach has been possible since the availability of 
high speed digital computers in the middle 1950s, but either was not 
attempted or was unsuccessfully attempted until the development reported 
on here . 

The earliest and only past effort in numerical-objective HBZ 
analysis known to this writer^, a statistical approach, was not success- 
ful, apparently because of an over specification of front-associated 
parameters and criteria for front location and intensity. It is a natural 
tendency to presume that the utilization of a multiple number of variables, 
such as pressure, wind, temperature, humidity, and their derivatives, would 
fare better than any one of them in objectively placing fronts. But, the 
diverse reactions of each of these parameters to the synoptic and mesoscale 
processes which affect baroc linicity , as modified by topography, compli- 
cates the standardization of locating HBZs by statistical formulae or by 
dynamical-statistical means. The more recent numerical-objective efforts 
by the United States Naval Weather Service (Hughes, 1967) and Air Weather 
Service (United States Air Force, 1965-1968), the Canadian Meteorological 
Service (Creswick, 1967), and British meteorologists (Kirk, 1965), have 
taken advantage of the research described in detail here and reported on 
earlier as preliminary accounts to the meteorological community (Renard 
and Clarke, 1965a, 1965b, 1966; Clarke and Renard, 1965, 1966a, 1966b; 
Renard, 1967a, 1967b, 1967c). 



Correspondence with personnel at the Third Weather Wing, United 
States Air Force, Offutt Air Force Base, Nebraska, and Travelers Research 
Center, Hartford, Connecticut; also, a Report Bibliography prepared for the 
writer by the Defense Documentation Center (1965). 
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Thus, the middle 1960s found numerical techniques and high-speed 
computers used to analyze, diagnose, and prognosticate the operationally- 
usable structure of the atmosphere, from near the surface to stratospheric 
levels, with one notable exception, namely frontal phenomena. Without 
philosophizing as to the inherent merits and deficiencies of frontal anal- 
ysis and prognosis to weather processes (Sawyer, 1964), their importance is 
evident immediately from the most casual perusal of daily weather maps from 
the United States (civilian and military) or any other nation. Such charts 
always include front locations for at least one level in the troposphere. 

An example, clearly illustrating the need for unique placement of 
fronts is shown in Fig. 1. This illustration displays the 0000 GMT 5 March 
1964 surface fronts analyzed by sixteen United States and foreign civilian 
and military weather centers. The non-uniformity of concept, definition, 
and model are manifested by fronts which cross at angles as much as 90 
degrees and differ in position by 300 miles or more when apparently indi- 
cating the same HBZ boundary. Such a situation is intolerable to a scien- 
tific discipline and, practically speaking, could be disasterous to a 
civilian or military operation which depended on such analyses for its 
success. What are the parametric criteria for placing each of the fronts? 

Are they concerned with temperature and/or temperature gradient, cloud and/or 
precipitation, vector wind and/or wind shift, pressure gradient and/or pres- 
sure tendency? Quantitative answers to these questions do not exist from 
any of the weather services whose frontal analyses are displayed in Fig. 1. 
The situation in 1970 differs very little from that in 1964 or 1960 (Jacobs- 
Haupt, 1963). 

In order to improve the matters just described, experiments were 
begun in 1964 to develop a model for numerically specifying hyperbar oc linic 
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zones, their boundaries and intensities, with the aim of objectively 
specifying these phenomena on a real-time operational basis. During the 
initial phases of the investigation, the decision was made to adopt a 
definition for the HBZ which is simple and acceptable to theoreticians 
and synopticians alike. Such a definition is that implied for the 
frontal zone in the American Meteorological Society *s Glossary of 
Meteorology (Brooks et al., 1959). In essence, the frontal zone is 
described as a three-dimensional zone or layer of large horizontal 
density gradient bounded by surfaces (i.e., frontal surfaces) across 
which this gradient is discontinuous. The intersection of a frontal 
surface with another identifiable surface or plane (e.g., a constant- 
pressure surface or the x, p plane) generates a line called the front. 

Such a definition agrees with that proposed by Taljaard et al. (1961) in 
their excellent appraisal of the state of frontal analysis in the early 
1960s. Theirs is still a very timely estimation. Since the concept 
stated in the Glossary of Meteorology clearly isolates density gradient, 
and thus temperature gradient, as the essential front- location parameter, 
the measure of frontal phenomena and baroclinicity are of the same nature. 
Therefore, the hyperbaroc linic zone, a bounded region of relatively large 
baroclinicity, is alternately referred to as a frontal zone and the HBZ 
boundaries as fronts. However, the term front is not limited only to the 
warm-side boundary of HBZs as conventionally analyzed in connection with 
frontal wave cyclones of the middle and high latitudes. Also, at this 
point, the discontinuity aspect of the frontal-zone definition should be 
interpreted as a mathematical specification not necessarily substantiated 
in the real atmosphere due to a limited data network and the scale of 



analysis . 
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Fig. 1. Surface fronts at 0000 GMT 5 March 1964, as trans- 
cribed from analyses of the following international and national 
civilian and military weather-analysis centers: 

National Meteorological Center, United States Weather Bureau; 

Central Analysis Office, Meteorological Branch, Department of 
Transport, Canada; Japanese Meteorological Agency; British 
Meteorological Office; Zentralamt, Deutscher Wetterdienst ; 

Icelandic Weather Bureau; United States Navy Weather Centrals at 
Alameda, Guam, and Pearl Harbor; United States Navy Weather Facilities 
at Argentia, Miami, Norfolk, Quonset Point, San Diego, Sangley Point, 
and Yokosuka. 



7 







8 



The subsequent chapters describe the development, characteristics, 
and applications of a numerical model for locating and measuring the 
intensity of real-atmosphere HBZs . Hemispheric, synoptic-scale case 
studies provide a limited four-dimensional (x, y, p, t) diagnosis of 
HBZs on a multi-seasonal basis. However, the reader may first wish to 
peruse Appendix I for a survey of the writer *s preliminary experiments 
in numerical frontal analysis. These first probes, somewhat air-mass 
oriented and based on wet bulb and equivalent potential temperatures, 
were not extended due to the difficulties of adequately computing and 
portraying moisture on a hemispheric scale as well as the artificiality 
of air-mass labeling. Relegating this section to an appendix signifies 
the terminal nature of these early experiments. Following these prelim- 
inary experiments, research was directed toward a more dynamic approach 
to numerical HBZ analysis, utilizing virtual potential temperature and 
its derivatives. The model proposed for general operational use and 
diagnostic testing is discussed at length in Chapter II. Chapter III 
illustrates the model by examples derived from hemispheric application 
in a framework of an existent operational analysis and forecast complex, 
namely the Fleet Numerical Weather Central, Monterey, California. A more 
refined scheme for specifying HBZs, featuring high resolution derived 
from vertical soundings, is presented next (Chapter IV) as a prototype 
for analysis on the subsynoptic scale. Chapters V-VII return to the 
operational model appropriate to present-day synoptic -sea le analysis 
(Chapter II) and focus on substudies concerned with the spatial and 
temporal consistency of numerical HBZ analysis and the relation of these 
analyses to cloud and vertical motion patterns. Finally, the treatise 
concludes with a resume* of the investigative results (Chapter VIII). 



CHAPTER II 



HYPERBAROCLINIC-ZONE ANALYSIS: 

DEVELOPMENT AND CHARACTERISTICS OF A NUMERICAL MODEL 

The Mode 1 

The preliminary experiments described in Appendix A, although 
valuable, merely suggested a continuation of the search for a more 
universal approach to objective hyperbaroc linic -zone analysis, that 
is, one void of intimate association with air-mass identification. 

The definition of the HBZ , basic to the analysis model whose 
development is outlined in this chapter, is described in Chapter I. 

Since the essential parameter is atmospheric baroc linicity , the 
numerical specification of the HBZ proceeds from the definition of 
the magnitude of total baroc linicity about a local point, which is 
given by Saucier (1955) as 

N = gjv In 0*1 . (1) 

N represents the number of unit solenoids per unit area in a plane 
determined by the near-vertical three dimensional pressure gradient 
and a normal to the isolines of virtual potential temperature (0*) 
on a constant pressure surface; g is gravity. Since V0* has mathe- 
matical properties similar to V In 0* while offering greater facility 
in handling of data, the former is used in the development that follows. 
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Symbolizing the unit vector normal to the line bounding the HBZ 
as n^ , a parameter for locating the line on a constant pressure surface 
may be defined by 



DG0*, = - V V0* -n, = 






5n 



2* 



( 2 ) 



Thus 



, DG0*^ is a directional second derivative of 0^, constrained to lie 



a long n^ . 

With reference to the identity in Eq. 2, the direction implied 

by n^ is unknown initially, as the position and orientation of the 

front is the quantity being sought. However, since the HBZ tends to 

— » 

follow an isotherm channel, a suitable approximation to n^ is given by 
the unit vector perpendicular to isotherms. Therefore, with a view 
toward quantitative applications, the HBZ-location parameter on a 
constant pressure surface is redefined as 



DG0* = - vlv0*| = - VG0*-nQ^ = 






0 * 



5n 



2 ^ 



(3) 



0^ 



where n^^ is a unit vector along V0*, and 

G0* = |V0*| = S0*/driQ^. 



(4) 



A One -dimensional Schematic 
Analysis of DG9* 

A one-dimensional representation of 0*, G0*, and DG0* which 
simulates the real atmosphere is an aid in describing the synoptic and 
thermo-physical properties of DG9*. Fig. 2 presents continuous 
schematic profiles of 0^ and its derivatives oriented along V0* on an 
isobaric surface. The profiles are slightly asymmetric about location 
0. Points 1 and 2 in Fig. 2_a are locations of relative maxima in the 
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Fig. 2. A continuous schematic profile of 0* and the G0* and 
DG0* profiles derived from it. 

a. 0*; theoretical profile with first-order discontinuity also 
shown in the vicinity of points 1 and 2 (dashed line). 

b. G0*5 derived from Fig. a. 

c. DG0*5 derived from Figs, a and b. 
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curvature of the 0* profile; theoretically, such points are first-order 
discontinuities in 0*. 

The expression G0‘*’ is a measure of the baroc linic ity (Fig. 2^), 
which is maximized at location 0 in the 0* profile. The inflection 
points (1, 2) in G0* are suitable boundaries of the zone containing 
relatively large magnitudes of baroc linic ity , i.e., the HBZ. A 
higher-order derivative of 0* is required to give a clear and quanti- 
tative indication of the boundaries of the HBZ. This is shown by the 
DG0* profile in Fig. 2c. The maximum and minimum values of DG0* 
indicate warm-side and cold-side boundaries of the HBZ, respectively, 
with the zero point representative of the location of maximum baro- 
c linic ity (Fig. 2b). Though magnitude of baroc linicity is not 
explicitly given by the front - location parameter, DG0*, a relatively 
large spatial change in V0* implies the existence of high values of 
baroc linic ity in the associated HBZ. Statistical and synoptic veri- 
fications of this relation are cited in this and later chapters. 

Physical significance of DG0* may be seen from the relations 
of geostrophic thermal wind, V^, and the isobaric virtual temperatures, 
T* and 0*. 



V I = 1.^1 = -S- ^ 9* 



fT* Sn, 






f 9n 



0 * 



and 






' hn 



0 * 



f 2 • 



( 5 ) 



( 6 ) 



Therefore, the shear component of the thermal vorticity, given by the 
shear of the thermal wind in Eq . 6, also exhibits a maximum in magnitude 



at each of the HBZ boundaries . 
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Analysis Characteristics of the Model 

In order to demonstrate the effects of numerical differen- 
tiation and interpolation, and scale on the location and intensity 
of HBZs , three hypothetical, but realistic, 0* profiles are presented 
(Fig. 3) and interpreted. The profiles are assumed to be oriented in 
the direction of V0*, here considered to be colinear with the I axis 
of the grid. Data are specified for whole -numbered grid points only 
(Table 1), consistent with numerically-analyzed observations. Although 
the distance between grid points, d, is somewhat arbitrary for the 
arguments set forth in this chapter, it will be of great practical 
value for what follows to consider d = 381 km^. 

The temperature scale in Fig. 3 is symbolic only, but as in 
the real atmosphere, the warm air is shown to be less baroclinic than 
the cold air in all three cases. The profiles are linear in the warm 
air but depart slightly from such in the cold air. The hyperbaroclinic 
zones, separating the coldest from the warmest air, are bounded by 
break points typical of nature *s version of first-order discontinuities 
on the synoptic scale. These maximum-curvature points in the modeled 
profiles are labeled as F^ and F^^ for boundaries on the warm- and cold- 
air sides of each of the HBZs, respectively. The profiles are labeled 
N (narrow), M (medium) and W (wide), descriptive of HBZ widths of one, 
two, and three mesh lengths, respectively. However, for all three pro- 
files, the HBZ A0* is 10 C. This results in HBZ gradients of 10 C 
(381 km)’^ for N, 5 C (381 km)”^ for M and 3.3 (381 km)"^ for W. 

^381 km is the mesh length at 60 deg lat in the FNWC grid 
representation of the Northern Hemisphere. 



r40 
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Fig. 3. Modeled profiles of 0* (solid lines), specified as N, M, and W, to identify 
the included narrow, medium, and wide HBZ widths, respectively, contained between F and F . 
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The aim of this chapter is to simulate, in a systematic but 
simple fashion, numerical differentiation and interpolation of tempera- 
ture data using computations of the type common to meteorological analysis. 
The modeled thermal data are considered numerically processed grid-point 
temperatures. Cases of definite physical-synoptic significance are sel- 
ected for simulation, as, for example, numerical processing in relation 
to intensity of frontal zones, mesh length, closely adjacent fronts, 
frontal slopes, and the use of layer vice level in the DG0* computations. 

The Simulation Experiments 



Frontal-zone Intensity 

Variation in HBZ width --First , the 0* data in Fig. 3 are used to 
generate G0* values by employing the Lagrangian 3- and 5-point numerical 
differentiation formulae (Jennings, 1964) for centered differences. For 
the 3-point computation at the central grid point, I^, the formulation is 



G0* = 
o 



2d 



( 7 ) 



while the 5-point G0* is given by the following expression: 



G0* = 
o 



8 [0^ - 0* J - 0^^. + 0* J 

o-t-1 o- 1 0+2 0-2 

12d 



( 8 ) 



G0^ grid-point data are then used to generate the DG0*s , again using 

both the 3-point and 5-point numerical differentiations. 

In order to show relative magnitudes of frontal parameters in 

and near the N, M, and W frontal zones. Table 1 lists the grid-point 

0^s (from Fig. 3) and the computed values of G0* and DG0^. It is to 

be noted that the magnitudes of G0* and DG0* at and between F and F, 

® w k 



Table 1. Grid point values of 0* (C), and computed G0* (C [ 100 km]‘ ) and DG0* (C[100 km]’ ) 
at intervals of II for the N, M, and W profiles in Fig. 3, both for the 3- and 5-point numerical dif- 
ferentiation schemes. Subsets of DG0* values indicate interpolated HBZ-boundary locations and their 
attendant DG0* maxima and minima (circled numbers). HBZ A0* = IOC. F (F, ) is the position of the 
warm (cold) boundary of the modeled HBZs . ^ 
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are similar to those found in frontal zones on the daily operational 
charts (Chapter III). 



directed is as follows. Assuming the modeled positions of the HBZ bound- 
aries to be at F and F, , where are the frontal positions as located by 
numerical differentiation and interpolation of the grid-point temperature 
data? To answer this question, Bessel*s central difference polynomial 
interpolation formula (Jennings, 1964) is utilized to interpolate the 
location of DG0* maxima and minima to the nearest .1 mesh length (i.e. 
about .3 deg lat for d = 381 km). These points are representative of the 
warm-side boundary (conventional front) and the cold-side boundary of the 
HBZ , respectively. 



where subscripts -1, 0, 1, and 2 indicate grid-point locations and s is 



The foremost question to which the profiles and tabular data are 



Setting DG0* = y, the formulation for DG0^ is 





+ - 2(y^-yo) + (y^-y,-^)]. 



the fractional distance between 1=0 and 1 = 1. The constants, Bj^, 



and B^ are functions of s in the form 



B 



1 



1 



and 




( 10 ) 



Consistent with the number of grid points used in the 5-point numerical 



differentiation formula (Eq. 8) the interpolation formula (Eq. 9) is 
truncated with the term which is cubic in s (third degree polynomial). 
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The numerically-interpolated locations of the HBZ boundaries 
and their attendant DG0^^ values are also shown in Table !• For 
example, in the case of the N profile the interpolated frontal-zone 
boundaries are interpolated to be at 7.4l and 5,4l (7.31 and 5.51) 
for the 3-point (5-point) computation with DG0^^ values of +0.21 and 
-0.17 (+0.32 and -0.27), respectively. By comparison, the modeled 
N-profile data indicate HBZ boundaries at 17 and 16 with DGO** values 
of +0.30 and -0.23, respectively. 

The following conclusions are pertinent to Table 1: 

1. The computed narrow (i.e. strong) frontal zone is broadened 
most relative to the modeled zone width, both for the 3- and the 5-point 
differentiation schemes. In fact, the weak frontal zone is fore- 
shortened, not broadened, in the 5-point case. 

2. The 3-point differentiation scheme broadens the HBZs more 
than the 5-point approximation. Specifically, the computed HBZ widths 
are broadened by the following percentages relative to the modeled zone 
widths ; 





N 


M 


w 


3-pt 


100 


15 


0 


5-pt 


80 


10 


-20 



3. The broadening effect is up to 20 percent greater on the 
cold-side boundary of the HBZ. It is to be noted that the cold-side 
boundary borders on the most baroclinic of the two air masses. 

4. The maximum grid-point magnitudes of 5-point G0* values 
exceed those of the 3-point scheme and may be greater than the values 
implied in the modeled HBZs of Fig. 3, as noted below for figures 



extracted from Table 1. 
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G0* 


(C[100 km]"^) 




HBZ 


Mode led 


Computed 








5-pt 


3-pt 


N 


2.63 


1.79 


1.58 


M 


1.31 


1.47 


1.31 


W 


.87 


1.00 


.92 



Hence, the 5-point G0* values are more desirable in order to yield 
frontal DG0*s which clearly outline the HBZs, 

5, The maximum DG0* is greater in magnitude than the minimum 
DG0* for a given HBZ , a consequence of greater baroclinicity in the 
colder of the two air masses. 

Further, note that if the modeled HBZ boundary in Fig. 3 is 
located between grid points the broadening and magnitude problems are 
increased. For example, in the case of the N profile, the broadening 
is quantitatively increased up to nine-tenths of a mesh length by 
placing the HBZ boundaries at 6.01 and 6.11 while maintaining the 0* 
value at 7.01. Thus, with a HBZ width of O.ll in the N case, the 
broadening percentage increases to 1700 percent for the 5-point dif- 
ferentiation. However, this circumstance is highly unlikely since a 
HBZ width < 40 km is unrealistic on the scale of present-day obser- 
vations and analyses. In case of the M and W frontal profiles the HBZs 
could be reduced to a zone lying between 5.91 and 6.11 and 4,91 and 
6.11, respectively, without a change in grid-point temperature data. 
This would result in a percentage broadening of 1000 percent and 100 
percent, respectively, for the M and W profiles in the 5-point case. 
However, nonmeteorological 0* profiles would result for the M and W 
cases if the grid-point temperature data remained unchanged. 
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Variation in HBZ AQ* - -Fronta 1-zone intensity (i.e. average A9* 
within the transition zone) is dependent on the change of 0* in the 
zone (i.e. A0*) as well as distance (An^^) through which A9* occurs. 

In order to demonstrate the effect of variations of V0* in the HBZ due 
to changes in A9*, the A9* is reduced to 5C for the N profile only of 
Fig. 3. This results in a modeled transition-zone gradient equivalent 
to that of the M profile in Fig. 3. Table 2 shows the 9*, G0*, and 
DG0* tabulations and broadening effects for the 3- and 5-point differ- 
entiations in the same format as Table 1. 

The broadening effect on the cold-side is greater than that on 
the warm-side boundary of the HBZ, an effect similar in kind to the N 
profile (A9* = IOC) even though the computed G0*s and DG0*s resemble the 
W case of Table 1. The overall percentage broadening is about the same 
as the N case in Fig. 3, namely 110 percent for the 3-point and 70 per- 
cent for the 5-point differentiations, again showing the desirability of 
the latter. In view of these effects, it appears that it is the frontal- 
zone width more than the magnitude of A9* which fosters significant 
computational broadening of the HBZ. That is, the real-atmosphere 
narrow transition zones are broadened relatively more than the wide 
zones by numerical processing, regardless of the magnitude of A9* across 
the HBZ. 

Both Tables 1 and 2 suggest the insignificance and noisiness of 
DG0* data outside the HBZ where values generally become < | .05|c(100 km) 
within one mesh length of HBZ boundaries. Also, the zero DG0* values lie 
half-way between the maximum and minimum values (to the nearest .1 d) for 
all 5-point cases in Tables 1 and 2. These facts provide guidance in the 
development of an operational DG0* product, as discussed in Chapter III. 
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Table 2. Grid point values of 0* (C), and computed G0* 

(C[100 km]"^) and DG0* (C[100 km]"^) at intervals of II for the N 
profile in Fig. 3, for the 3- and 5-point numerical differentiation 
schemes. Subsets of DG0* values indicate interpolated HBZ-boundary 
locations and their attendant DG0** maxima and minima (circled numbers). 
HBZ A0* = 5C. F (F, ) is the position of the warm (cold) boundary of 
the modeled HBZsY ^ 
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Variation in 
Mesh Length 

In order to demonstrate the increase of resolution in locating 
fronts with a decrease in mesh length, the N and M profiles of Fig. 3 
were considered in connection with a doubling of the number of grid 
points. Linear interpolation was used to obtain the 0* values at the 
positions intermediate to the whole numbered grid points in Fig. 3 (i.e., 
at ^ 2 etc.) with slight exception in the M case. Here, a 

deviation from linearity was employed to show maximum baroc linicity at 
the geometric center of the HBZ . The 5-point differentiation scheme 
and A0* = IOC are common to the two cases computed. 

N Case --Table 3 gives the pertinent information on this case. 

The .positioning of the HBZ boundaries is precise, hence the broadening 
effect is zero compared to 80 percent broadening with data at whole- 
numbered grid points only. The maximum G0* value has increased by 
almost 70 percent compared to the N case in Table 1 with G0* exceeding 
the model value of 2.63 C (100 km) the associated DG0*s have tripled, 
a desirable result for quantitative specification of the HBZ boundaries. 

M Case --Table 3 also gives the details on the M case. Much like 
the N case the positioning is excellent, with the HBZ foreshortened by 
two-tenths mesh length (10 percent) compared tothemodeled width, mostly 
due to conditions on the cold-side boundary. Again, magnitudes of 
and DG0* have increased compared to their counterparts in Table 1. 

The obvious conclusion is that a mesh-length reduction resolves 
the computed position and intensity of the HBZ. Any nominal fore- 
shortening due to numerical processing is not only tolerable but desir- 
able due to the smoothing inherent in generating the data at grid points 
from the reported observations . 
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Table 3. Grid point values of 0* (C), and computed GG'* 

(C[100 km]*i) and DG9* (C[100 km]*^) at intervals of .51 for the 
N and M profiles in Fig. 3, for the 5-point numerical differentiation 
scheme only. Subsets of DG0* values indicate interpolated HBZ-boundary 
locations and their attendant DGG^ maxima and minima (circled numbers). 
HBZ A0* = IOC. F (F ) is the position of the warm (cold) boundary of 
the modeled HBZs ^ 
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Closely Adiacent Fronts 

A typical occurrence is that which may be described as closely 
adjacent fronts. Operationally, the fronts are referred to as primary 
and secondary; the former designates the most equatorward of the two 
fronts. Usually both frontal zones are of the cold-front variety and 
the HBZ V0* associated with the secondary is generally weaker than that 
of the primary. 

Two examples are given to illustrate this case. (1) A combin- 
ation of N and M profiles with an intervening subbaroclinic zone or 
air mass of one-mesh length dimension; appreciable baroc linicity occurs 
in the air separating the N and M transition zones. (2) Same as (1) 
except an intervening air mass of two mesh lengths. 

Table 4 illustrates the multiple-front problem with values of 
G0* and DG0*. 

One-mesh air mass --In the one-mesh air-mass case the broadening 
effects noted with N and M profiles (Table 1) are reduced. In other 
words, the N or primary frontal zone shows 40 percent vice 80 percent 
broadening and the k-side is more accurately positioned than the w-side, 
compared to the uncomplicated profiles in Fig. 3. The M or secondary 
frontal zone is foreshortened by 40 percent (vice 10 percent broadening 
in Table 1) with the k-side again placed most accurately. The secondary, 
especially its w-boundary, behaves as a very weak front here. 

When viewing the intensities, the w-boundary of the primary is 
similar to profile N in Table 1, while the k-boundary is reduced in 
magnitude. The w boundary of the secondary is close to the magnitude of 
the noise level criterion described earlier, hence the secondary is 
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Table 4. Grid point values of 0*^*^ (C), and computed G0* 

(C[100 km]’^) and DG0* (C [ 100 km]"^) at intervals of II for the 

multiple-front case with air-mass dimensions of one and two mesh 

lengths. The primary (secondary) front is represented by the N (M) 

HBZ profile in Fig. 3; 5-point differentiation scheme only. Subsets 

of DG0* values indicate interpolated HBZ-boundary locations and their 

attendant DG0* maxima and minima (circled numbers). HBZ A0* = IOC. 

F (F ) is the position of the warm (cold) boundary of the modeled HBZs . 
w k 





1-Mesh 


Air Mass 






2 -Mesh 


Air Mass 




I 


gic 


G0* 


DG0* 


I 


e* 


G0* 


DG0* 


13 


32.55 






13 


32.55 






12 


32.5 






12 


32.5 






11 


32. 


.14 




11 


32. 


. 14 




10 


31.5 


.13 




10 


31.5 


.13 




9 


31. 


.13 


-.04 


9 


31. 


.13 


-.04 


8 


30.5 


.08 


+ .21 


8 


30.5 


.08 


+ .21 


7 


30. 


1.56 


+2.8F 

w 


7 


30. 


1.56 


+ .30F 

w 


6 


20. 


1.72 


-. 17F, 
k 


6 


20. 


1.79 


-.26F, 

k 


5 


18. 


.74 


-.03F 

w 


5 


18. 


.21 


-.17 


4 


13. 


1.46 


+ .04 


4 


16.5 


.84 


+ .24F 

w 


3 


8. 


.79 


- . 2 IF, 
k 


3 


11.5 


1.46 


+ .01 


2 


7. 


.15 


-.08 


2 


6.5 


.94 


-.18F, 

k 


1 


6.25 


.13 


-.01 


1 


4.5 


.33 


-.10 


0 


6. 


.02 




0 


3.5 


.21 


-.01 


-1 


6. 


.01 




-1 


2.75 


.20 




-2 


6. 






-2 


2. 


.16 




8 






+ .21 


8 






+ .21 


7.3 








7.3 






<£B> 


7 






+ .28F 

w 


7 






+ .30F 

w 


6 






-.17Fj^ 


6 






-.26F 


5.9 








5.7 






dD 


5 






-.03F 

w 


5 






-.17 


4.2 






(JTg) 


4 








4 






+ .04 


3 






+ .01 


3 






ci>-k 


2 









26 



poorly resolved. However, the secondary’s k-boundary is well defined 
in magnitude. 

It should be noted that a one-mesh air-mass distance is not 
very realistic. Further, as in the simple structures in Fig. 3, the 
numerical processing effect causes all HBZ widths to approach a medial 
value. The relatively narrow frontal zones broaden, the wide ones tend 
to foreshorten. 

Two-mesh air mass --When the air-mass zone is widened to two 
mesh lengths both the primary and secondary zones behave as if the 
other were not present. In the two mesh case, the primary HBZ is 
broadened by 60 percent vice 80 percent while magnitudes are approxi- 
mately the same as those of the N profile in Table 1. The secondary 
front is precisely located and magnitudes agree closely with the simple 
M profile in Table 1. 

It is obvious from this experiment that the position and intensity 
of each of the HBZs become increasingly independent as the subbaroclinic 
zone separating them widens . 

Vertical Slope of 
HBZ Boundaries 

What is the relationship of the slope of the HBZ boundaries and 
the vertical variation of HBZ widths as caused by numerical processing? 
There are two cases of interest; transition -zone baroclinicity weakens 
with increasing height; and transition-zone baroclinicity intensifies 
with increasing height. These are applied to the simple profiles (N, M, 

W) shown in Fig. 3. The frontal-zone intensity is a function of the 
width only. 
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HBZ weakens with increasing height --Fig. 4a illustrates the 
situation for a modeled warm-boundary slope of approximately 1/240, 
with frontal-zone characteristics changing from N to M to W from 1000 
to 850 to 700 mb, respectively. This results in an associated cold- 
boundary slope of 1/480. The 1000 to 700 mb warm-boundary slope, 
determined from numerically-computed data, is 1/312, which is a 23 
percent reduction from the modeled slope. At the cold boundary the 
average slope increases relative to the hypothetical slope by 26 per- 
cent (from 1/480 to 1/384) in the layer 1000 to 700 mb. 

The computed transition-zone widths, in comparison with the 
modeled values, range from greater (i.e. broadening) at 1000 mb to 
less (i.e. foreshortening) at 700 mb. The overall effect is that of 
reducing the range of HBZ widths relative to the model values. Boundary 
values of DG0* decrease with height consistent with the vertical gradient 
of transition-zone baroc linicity . 

HBZ intensifies with increasing height --Fig. 4^ illustrates the 
case of the frontal intensity changing from W to M to N with increasing 
height in the layer 1000 to 700 mb. Again, 1/240 is taken as the slope 
of the warm-air edge of the HBZ; the computed value is 1/168. The 
modeled cold-side boundary of the HBZ is constrained to be vertical with 
the computed slopes ranging from 1/120 in the low layer to 1/72 in the 
upper layer. Thus, both the broadening effect and magnitude of DG0^ 
maxima increase with increasing height. 

Use of Layer-mean Temperature, 
in the DG0^ Analysis 

Low-tropos phere mean temperature (i.e. is obaric- layer thickness) 
fields enjoy widespread use for locating fronts operationally. Here, 
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Fig. 4. Numerically-computed HBZ boundaries (dashed lines) in 
relation to the modeled boundaries (solid lines); HBZ A9* = 10 C. 

a. Slope of modeled warm-side boundary: 1/240; width of HBZ 

increases with height 

b. Slope as in a; width of HBZ decreases with height 
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typical 1000/700 mb mean- layer temperature (0*) profiles are generated 
from information in Fig. 3 and reasonable vertical slopes. Only the 
5-point derivative scheme is used. The following cases are considered 
for HBZ slopes of 1/120 and 1/240: (1) 1000 and 700 mb --N-type HBZs 

both levels; (2) 1000 mb --N-type HBZ, 700 mb --M-type HBZ; and (3) 

1000 mb --M-type HBZ, 700 mb --N-type HBZ. 

0*, at each grid point, is obtained as an arithmetic average of 
the 1000- and 700-mb information where the appropriate 700-mb profile 
(N or M) is displaced one (two) mesh length(s) toward the cold air for 
the 1/120 (1/240) slope case, respectively. In addition, 3C to 8C are 
added to the 0* values at 700 mb in order to simulate the normal increase 
of 0* with height, with the higher (lower) temperatures used in the warm 
(cold) air mass. The addition is made in such a way as to preserve the 
transition zone V0* at each isobaric level. Such a manipulation of the 
data agrees with the observable result that the cold air is both more 
unstable and baroclinic than the warm air in the lower part of the 
troposphere. Table 5 shows 0*=*^, G0^, and interpolated HBZ boundaries 

similar to previous sections. 

N-type, both at 1000 and 700 mb; slope: 1/120 --This case simu- 

lates a relatively strong cold front with steep slope. The numerical 
positioning of the HBZ boundary is exact, resembling the M profile, 
single level, while the magnitudes of DG0* are intermediate to those 
for N and M at a single level (Table 1). The HBZ width for 0* is the 
distance from the 1000 mb warm-side to the 700 mb cold-side HBZ boundary, 
thus essentially doubling the hypothetical zone width compared to either 



1000 or 700 mb. 
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Table 5. Gr^d point values of 1000/700 mb 0*, and computed G0* 
(C[100 km]‘l) and DG0* (C[100 km]"2) at intervals of ll for N, M, and W 
profile data derived from Fig. 3 and frontal slopes__of 1/120 and 1/240; 
5-point differentiation scheme only. Subsets of DG0* values indicate 
interpolated HBZ-boundary locations and their attendant DG0* maxima and 
minima (circled numbers). HBZ A0* = IOC. position of the 

warm (cold) boundary of the modeled HBZs . 
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The nature of the vertical slope is such as to give the effect 

of an M case in both frontal intensity and width since the superposition 

of F, 1000 and F 700 tends to obliterate the identity of each of these 
k w 

HBZ boundaries when using the mean temperature field. 

N-type both at 1000 and 700 mb; slope: 1/240 --Pos itionlng is 

poor, with the numerical 1000 mb front displaced by one-third mesh length 
toward the warm air and the computed 700 mb cold-side HBZ boundary ex- 
tended back into the cold air by eight-tenths mesh length, a 36 percent 
increase over the modeled HBZ width. Magnitudes of DG0* are reduced to 
those of the W category, but the broadening is typical of the N category 
(Table 1), partly due to the greater geometric separation between the 
1000 mb warm-side and 700 mb cold-side of the HBZ compared to the 1/120 
slope case. 

For further slope reduction the minor DG0* minima and maxima, only 

evident as noise level values at 16 and 15 in Table 5, become stronger, 

until two separate and significant HBZs , F 1000 to F, 1000 and F 700 to 
^ w k w 

F^700, result. Limiting such occurrences is the fact that in the real 
atmosphere, the N-type fronts have large slope (i.e. > 1/120). 

1000 mb--N-type; 700 mb--M-type; slope: 1/ 120 --Pos ition and 
intensity of the warm-side boundary of the HBZ behaves as the M type but 
the cold-side boundary, F 700, appears as a doublet with small magnitudes 

K 

in each case. The cold boundary complexity results from the superposition 

of a strong F, 1000 and a moderate F 700 front. 

^ k w 

1000 mb -- N-type ; 700 mb--M-type; slope: l/240 --In this case, 

due to increased separation of the HBZs at each level, the 1000 mb and 
700 mb frontal zones manifest themselves as separate entities. With 
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further reduction in slope the F, 1000 and F 700 boundaries become 
increasingly distinct. In this case, the F^700 value of +0.03 may 
be disregarded as noise. 

1000 mb --M-type; 700 mb --N-type; slope: l/l20 ""Some rather 

unexpected results show up for this case. The warm boundary is fore- 
shortened, whereas a broadening effect occurs at the cold boundary. 

Magnitudes are typical of the M case (Table 1). Apparently, the F 700 

w 

position is the cause of the peculiar positioning of the 1000 mb front 
giving the properties of the W profile to the warm boundary, while N 
properties are noted at the rear boundary where F 1000 and F 700 are 

tC K. 

superimposed . 

1000 mb ""M-type; 700 mb ~~N-type; slope: 1/240 - -Comment is 

similar to that immediately above with the added feature of a more 
intense cold-side HBZ boundary compared to the warm-side. 

Summary of the 1000/700 mb layer simulation --Use of 9* vice 0* 
appears to offer advantages for the front-location model. The broaden- 
ing effect at the HBZ warm-air edge is reduced, hence this entity is more 
accurately located by numerical means. Although DG0* values are less than 
DG9* for similar transition zones, the magnitudes remain well above noise- 
level values. With a view toward operational analysis of DG0*, it is not 
too serious that the cold-side boundary locations are complex or poorly 
located due to the smearing inherent in the thickness model since, by 
convention, only the warm edges of frontal zones are portrayed on the 
surface or near-surface isobaric charts. 

The occurrence of fronts intensifying with height tends to produce 



foreshortening (Table 5) of the HBZ, which helps to balance the broadening 
effect due to sparse data. 
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General Conclusions 

The general conclusions concerning the accuracy of numerical 
HBZ analysis as derived from experimentation with the profiles in Fig. 

3 , are as follows : 

1. Use of the Lagrangian 5-point vice the 3-point numerical 
differentiation formula results in the most accurate placement and 
intensity of the HBZ. 

2. There is a tendency toward uniform HBZ widths in the 
numerical analysis of DG0*, that is, the narrower zones are widened, 
the wider zones are narrowed. Stated another way, broadening of the 
HBZ is more of a function of scale than of degree of baroc linicity in 
the zone. 

3. The most intense or sharply-focused HBZ boundary is the 
warm-side (i.e., the conventional front), just as in the actual atmos- 
phere, a result of the relatively smaller amount of baroc linicity in the 
warm compared to the cold air. 

4. Although there is a definite tendency toward an unreal 
uniformity of computed transition-zone widths, there is difficulty in 
devising an empirical correction as a function of the magnitude of the 
DG0* maximum since, by itself, frontal intensity (V0*) is not a good 
indicator of the numerically-derived width. Rather, it is the true 
frontal-zone width which is a good indicator of the broadening effect, 
but true width is not known. However, it is apparent that the frontal 
zones with relatively large thermal gradients are most likely to be 
broadened but the weaker-gradient zones may be broadened or fore- 
shortened due to numerical processing. 
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5. Overall, the 0* model is more desirable than the 0* model 
since the effect of using a mean temperature biases the front - location 
properties toward those of the wide and weak-zone variety and this 
reduces the broadening effect. However, when using 0*, the cold boundary 
of the lower-level frontal zone cannot be located. This is also true, 

to a large degree, for the cold boundary of frontal zones at the upper 
level of the layer. 

6. Closely adjacent frontal zones are likely to be misplaced. 

In that case, foreshortening and weakening of the conventional boundary 
of the secondary is a likely effect. 

7. Reduction of mesh size aids in resolving the true location 
and intensity of HBZs . Any tendency toward foreshortening of the HBZ , 
as mesh size is reduced, is also regarded as favorable considering the 
smoothing inherent in data processing and the relative sparsity of data. 

8. The vertical variation of HBZ intensity, in terms of the 
zone*s width, contributes toward uniformity of widths and vertical slopes 
of the HBZs. 

It should be emphasized that the general analysis characteristics 
of the DG0* model presented in this chapter were reduced from a rela- 
tively simple set of temperature profiles and a controlled simulation 
experiment. The real atmosphere presents considerably greater com- 
plexity in temperature distributions with a wide range of concurrent 
conditions of frontal-zone intensity, width, slope, and contiguity. 

Thus, it is logical to consider next a two-dimensional application of 
the DG0* numerical analysis scheme to a hemispheric portion of the real 



atmosphere . 



CHAPTER III 



AN OPERATIONAL APPLICATION OF THE 
HYPERBAROCLINIC-ZONE ANALYSIS MODEL 

Chapter II summarized the development of a numerical analysis 
model which has potential use in identifying, locating, and diagnosing 
HBZs in real time. Insight was also given on the effects of scale and 
numerical differentiation and interpolation on the computed values of 
baroclinicity and zones of hyperbaroc linicity . To this point in the 
presentation, neither operational merits and deficiencies of the DG0* 
analysis nor the effects of numerical data processing on the scheme 
have been considered. Thus, it became desirable to program the DG0* 
analysis scheme for use in a computer-oriented weather central with 
international responsibility. Such an environment, and the one most 
accessible to the physical location of the subject research, is that of 
the United States Navy’s Fleet Numerical Weather Central, Monterey, 
California (FNWC) . There follows a discussion of the FNWC automatic 
data processing program and examples of two-dimensional DG0* analyses 
as performed in that framework. Since field use is directed toward the 
analysis of the conventional front, numerical-objective and manual- 
subjective versions of frontal analyses are discussed in this and later 
chapters . 
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Data-Processing Framework for the Operational DG0^ Analysis 
Germane to the research on hyperbaroclinic zones is the FNWC data 



processing framework in which the DG0* analysis is performed. Roll, et 
al, (1963, 1965) have developed a so-called mass structure model for 
numerically processing meteorological data in a physically consistent 
manner. The model is used currently by FNWC (Hughes, 1967). The fea- 
tures of the model, which relate directly to a comprehension of the 
peculiarities of the analysis of baroclinicity , are described below. 



analyzed by FNWC from a Northern Hemisphere synoptic- time coverage of 
more than 4000 surface and 500 radiosonde observation stations. These 
data are processed at 12-hour intervals to produce objective analyses of 
sea-level pressure and height and temperature of all mandatory pressure 
levels up to 200 mb, as well as dew-point depressions at terrain level, 
850, 700, and 500 mb. Except for terrain- level parameters, the analyses 
are produced on a square grid, 63 x 63, wherein the equator is an 
inscribed circle. 

The mesh length is 381 km at 60 deg lat on a Northern Hemisphere 
polar stereographic projection. At other latitudes, $, the geometric 



The mass -structure processing of input data is accomplished in 
the following manner. A mean static stability for each of five layers 
(1000 - 775, 775 - 600, 600 - 400, 400 - 300, and 300 - 200 mb), at each 
radiosonde station, is determined as a function of the heights of the 
mandatory pressure levels and sea -level pressure. Objective analyses 



The data used as input to the DG0* program are those objectively 




of the layer stabilities, the 1000-mb height, and the thickness of the 
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layer 1000 to 500 mb yields grid point values of these key parameters 
at each of the 3969 grid points. By a complex matrix solution, grid 
point values of temperature and height, at each of the mandatory pres- 
sure surfaces, are generated from the fields of the key parameters. The 
temperatures so produced are linear in p at each grid point (p = 

pressure; k = R,/c ; R, = the gas constant for dry air; c = the specific 
dp d p 

heat at constant pressure) and closely approximate virtual temperatures. 
The objectively-analyzed fields of the state parameters, pressure (or 
height), temperature, and moisture, are then available to compute other 
meteorological quantities, such as those selected for numerical frontal 
analysis . 

For computing derivatives, the grid-point values are approximated 
by employing the Lagrangian 5-point differentiation formula for centered 
differences along each of the linear grid axes, I, J, with the two com- 
ponents combined vectorially. See Eq . 8, Chapter II. As a rule, the 
first derivatives were processed with a low-pass filter, having a cut- 
off wavelength of four mesh lengths (approximately 850 nautical miles) 
before further differentiation. 

An Example of Two-Dimensional Analyses of DG0* 

Fig. 5 indicates the typical mathematical relationships 
existing between FNWC's objectively-analyzed isobaric fields of 0* and 
the patterns of G0* and DG0* derived from it for a 5 x 8 mesh- length 
area (d = 381 km). The illustrations may be regarded as simulations of 
middle- latitude winter conditions in the low troposphere. 

The plan views contain two separate zones of concentrated thermal 
gradients (i.e., HBZs), the bounding lines of which are shown in each of 
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the Figs. 5a_, and c. The field of 0* in Fig. 5a is used to generate 
the isolines of G0* in Fig. 5_b. The HBZ associated with the coldest 
temperatures is most baroclinic with G0* exceeding 2.0 C (100 km) ^ 
along the zone^s centrum. Fig. 5£ shows the DG0* field, exhibiting 
elliptically shaped isolines. 

The locus of points along which V0^^ changes most rapidly in the 
direction of n . defines the axes of maximal and minimal DG0*. These 
are the warm- and cold-side boundaries of the HBZ, respectively. Both 
boundaries tend to parallel the isentropes, especially the former (Fig. 

5^) . It is to be noted that the axes of maximum and minimum V0* (from 
Fig. 5b) nearly coincide with the zero isolines in Fig. 5£; these lines 
represent the centrums of the hyper- and sub-baroc linic regions, res- 
pectively. The imperfect coincidence of G0* axes and the zero isolines 
of DG0* is a result of the numerically-computed finite difference 
approximations to the analytic expressions. Fig. 5d illustrates pro- 
files of 0*, G0*, and DG0^^ along line N in Figs. 5£--£. The similarity 
to the schematic profiles in Fig. 2 is obvious. 

Operational Hemispheric HBZ Analysis 

Following closely on the development of the HBZ analysis model, 
FNWC has been producing, twice daily, operational hemispheric 1000 mb 
numerical frontal analyses for the guidance of Fleet Weather Centrals 
and Facilities in analyzing the important synoptic-scale low-troposphere 
fronts. Analysis criteria, degree of smoothing, and other presentation 
features have evolved during the period 1965 to date. 

In view of the experiments in Chapter II, a mean virtual potential 
temperature for the layer 1000/700 mb, 0*, was selected to generate the 
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Fig. 5. Typical plan views and profiles of 0*, G0*, and DG0*, 

a. 0* isolines, at intervals of 2.5 K. Positive and negative 
axes of DG0* superimposed, the latter as dashed lines, the 
former as dashed lines with alternating circles and diamonds. 

b. G0* isolines, at intervals of 0.2 C (100 km)"^. Axes in the 
G0* field shown as dotted lines. DG0* axes as in a, 

c. DG0* isolines, at intervals of .05 C (100 km)"^, with axes 
as in a; G0* axes from b, 

d. Profiles of 0*, G0*, and DG0* along line N in £, by and £. 
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front- location parameter appropriate to the 1000 mb (or surface) chart, 
DG0^. For this purpose, 0* is derived from the layer thickness and 
assigned to the logarithmic mean isobaric level, 836 mb. The next 
subsection focuses attention on the merits and deficiencies of the 
operational DG0* analysis with illustrations for a typical winter day, 

A statistical evaluation of numerical and manual frontal analyses serves 
to highlight differences and similarities of the two modes of analysis. 

Objective Frontal Analysis for 0000 GMT 6 March 1966 

The radiosonde observation time, 0000 GMT 6 March 1966, was 
selected to illustrate the features of objective frontal analysis. This 
particular synoptic time was not chosen because of the excellence of the 
numerical frontal analysis, but rather because of availability of many of 
the charts, fields, etc. (Figs. 6 to 11) necessary to demonstrate the 
present status of the analysis using the DG0* scheme. 

Fig. 6 shows a typical numerically-produced analysis of the 
frontal parameter on a polar stereographic projection of the Northern 
Hemisphere. Because layer temperatures are used and since the operational 

application is to (conventional) frontal analysis, only isolines of 

— -2 -2 
positive DG0* are shown in Fig. 6, starting with 5 x 10 C (100 km) 

Values of DG0* lower than 5 units fall in the noise range. The axes of 
relative maxima in the field of the elliptically-shaped isolines repre- 
sent the positions of the warm-side boundaries of the 1000 mb HBZs , hence 
they may be used by the field meteorologist as a first guess to locate the 
significant surface fronts . As an example, the numerically-analyzed 1000 
mb fronts in the middle and eastern Atlantic area are indicated by the 
superimposed dashed lines which trace out the DG0* axes. 
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Fig. 6. DG0* (1000/700 mb) analysis for 0000 GMT 6 March 1966. 
Isolines at interval of 20 x 10 "^ c (100 km)”^ , starting with +5 units. 
DG0* axes indicated by dashed lines for a section of the Atlantic area, 
NMC fronts are shown with conventional symbols. Alphabetical desig- 
nators A, B, C, D refer to Fig. 16. 
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What constitutes a significant front or HBZ boundary may be 
determined by the spatial extent of a DG0* axis and the attendant mag- 
nitudes, as well as by the location, general orientation, and shape of 
this field. The nearly meridional eastern United States-North Atlantic 
and nearly zonal western Pacific DG0* axes are examples of significant 
fronts in Fig. 6. Extensive zones with DG0* in excess of 25 units 
signify intense fronts; these are HBZs containing large temperature 
gradients. Note the northern Pacif ic--Siberia and western Pacific-- 
China fronts with DG0* exceeding 25 units over a large area; maximum 
values are greater than 65 units. 

The orientations, but not necessarily the positions, of axes 
of DG0* in Fig, 6 appear reasonable relative to the fronts usually seen 
on the daily synoptic analyses emanating from the National Meteorologi- 
cal Center, Suitland, Maryland (NMC) . However, air-mass designations 
(now mostly extinct, an}rway) and movement (warm, cold, etc.) are not 
indicated explicitly in the DG0* field; also, stage of development, as 
the occluded stage, is difficult to interpret from the numerical analysis. 

The problem of locating occlusions in the numerical frontal 
analysis is manifest over the eastern North America- -wes tern Atlantic 
region and the central-western Pacific. Minor bulges in the DG0* iso- 
lines may be associated with occluded structure but a unique empirical 
solution for this aspect of the analysis is difficult to program. The 
occlusion problem relates to the direction in which the virtual potential 
temperature derivatives are computed (Eqs . [2] and [3]). A typical view 
of NMC's subjectively-analyzed surface occlusions in relation to the 
1000/500 mb thickness lines serves to illustrate the problem well (Fig. 7). 
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The occluded fronts over the Pacific Ocean west of 180 deg and over the 
western and eastern United States intersect thickness lines (i.e., mean 
isotherms) at angles which, in many instances, are as much as 90 deg. 

Such frontal orientations are not likely to be derived by computations 
of DG0* from Eq . 3, as may be noted in Fig. 6. With reference to Eqs . 

(2) and (3), if the front - location derivatives could be calculated in 
the n. rather than in the n direction the problem associated with the 
numerical analysis of occlusions would not exist. A later section of 
this chapter discusses a modification to the DG0* model which somewhat 
alleviates the problem. 

The multiplicity and orientation of fronts implied in Fig. 6 has 
more similarity to the frontal-contour analyses performed by the Canadian 
Meteorological Service (Anderson et. al., 1955) than by NMC . Lacking an 
absolute standard, analyses from both weather services have been taken as 
the best available, especially in the dense-data North American and North 
Atlantic areas, and, as such, have been used in statistical comparisons 
with the objective frontal positions. A later section of this chapter 
presents results from the coevaluation of NMC and objective fronts. 

Numerically-analyzed fronts tend to be stronger over land than 
ocean; fronts near the land-ocean interface of continents (especially 
eastern continents) are pronounced. Part of the reason for these 
analysis features is due to the sparsity of data over oceans and areas 
of rugged terrain. However, it is also true that the relatively large 
intensity over land is due to the rapidly-varying height and nature of 
the underlying surface. Use of 0* vice 0* has not eliminated topographic 
effects on the objectively-derived frontal patterns, as indicated by the 
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analysis in the vicinity of the Rockies, Sierra Madres , Greenland, and 
Himalayan areas in Fig, 6. Studies of monthly-mean DG0* axes (Chapter 
VIII) and attempts to remove orographic influence and interface effects 
are germane to this problem (Reed, 1961). 

As indicated in Chapter II, limiting the DG0* plot to 5 units and 
greater screens out the noise frequently associated with the very low DG0* 
values in the subbaroc linic regions. Such a scheme also avoids dis- 
playing those insignificant zones of hyperbaroclinicity where DG0* values 
are < +5 units . 

The occurrence of a double frontal structure in NMC's analysis 
for the western Pacific but only a single HBZ boundary in the numerical 
analysis and the tendency for DG0* positions to occur on the warm-side 
of NMC positions, as shown both in the west and east Atlantic (Fig. 6) 
are manifestations of broadening effects and the close proximity of HBZs . 
Such problems were discussed in connection with the simulation experi- 
ments in Chapter II. Grossness of mesh length, peculiarities of the data 
processing, and the differentiation and interpolation schemes used 
play a part in explaining the differences in positions of manually- 
analyzed and numerically-analyzed fronts. 

Fig. 8 illustrates FNWC’s operational version of the DG0* analysis 
for 0000 GMT 6 March 1966. The operational chart is a function of the 
limited amount of data received by observation time plus four hours . This 
fact and the retention of a few small isolated DG0* zones in the opera- 
tional version (areas A and B) account for the minor differences existing 
between Figs. 6 and 8. 

Figs. 9 and 10 show the operational 1000 mb numerically-analyzed 
fronts superimposed on the sea- level pressure field and the 850 mb 
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Fig. 8. Operational DG0^ (1000/700 mb) analysis for 0000 GMT 
6 March 1966. Isolines at interval of 5 x 10*^ C (100 km)“2 , starting 
with the -♦‘5 line. DG0* axes (dashed lines) are manually positioned for 
portions of the Atlantic area, similar to Fig. 6. A and B indicate 
insignificant DG0* features. 
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Fig. 9. FNWC sea- level pressure field, 0000 GMT 6 March_1966. 
Isobar interval; 4 mb. 1000 mb numerica lly- analyzed fronts (DG0* axes) 
from Fig. 8 (heavy solid lines) with movement characteristics shown for 
the west and east Atlantic Ocean fronts by conventional symbols. 
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Fig. 10. FNWC 850 mb isotherm and contour field, 0000 GMT 
6 March 1966. Isotherms (thin solid lines) at 5 C interval, contours 
(thin dashed lines) at 60 met interval. 1000 mb numerically-analyzed 
fronts (DG0* axes) from Fig. 8 (heavy solid lines) with movement 
characteristics shown for the west and east Atlantic Ocean fronts by 
conventional symbols. 
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isotherm-contour chart, respectively. The symbols designating warm, 
cold, and stationary are derived from geostrophic winds at sea-level and 
thermal conditions at 850 mb. Most of the DG0* axes are located in or 
are close to troughs in the pressure field, giving evidence of internal 
consistency in the mass structure model used to process the input data 
at FNWC« The 850 mb isotherm pattern agrees quite well with the position 
of the objective fronts. There is a distinct tendency for the isotherms 
to parallel the HBZ boundary except near pressure centers, similar to 
manually-analyzed fronts (Fig. 7). 

Objective Frontal Prognosis for 1200 GMT 7 March 1966 

The numerical model described here may be used equally well for 
objectively locating fronts in the prognostic thermal fields. Fig. 11 
shows the operational 36-hour 1000 mb prognosis of DG0* and the verifying 
chart for 1200 GMT 7 March 1966. The FNWC prognostic model used to gen- 
erate the isotherms from which DG0* is computed is a modified barotropic 
scheme (Hughes, 1967). Such a model may be expected to smooth or elim- 
inate the minor thermal waves but leave the more intense and extensive 
synoptic -sea le features intact. It is for this reason that the leve 1 
(DG0*) vice the layer (DG0*) HBZ analysis is applied to the prognostic 
thermal fields. It is emphasized that the prognostic frontal positions 
result from an ana lysis of prognostic thermal fields. 

Comparison of Subjective and Numerical Frontal Analyses 
1000 mb Analyses 

Comparative statistics on the objectively- and subjectively- 
analyzed frontal positions serve to quantify the observations on DG0* 
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Fig. 11. a. Operational 36-hour prognosis of 1000 mb DG0* for 1200 GMT 7 March 1966. 
Interval; 5 x 10'2 C (100 km)'^, starting with the+5 isoline, 
b. Verifying 1000 mb DG0* analysis for 1200 GMT 7 March 1966. Interval; 

5 X 10"2 C (100 km) "2, starting with the+5 isoline. 
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made earlier in the chapter as well as give direction to the search for 
modifications to the present program. It should be emphasized that an 
absolute standard for locating HBZs does not exist, therefore it is 
necessary to compare manual frontal analyses of high quality with the 
numerical product. The NMC analyses were chosen for this purpose. 

The most extensive evaluation was made for a late summer period 
of sixteen observations times from 23 August to 6 September 1965 (Table 6). 
Over 250,000 nautical miles of numerically-located 1000 mb fronts were 
compared with the analogous NMC surface fronts. The more obvious con- 
clusions from this study are as follows: 

1. The two independent frontal analyses have the closest relation 
over dense-data land and ocean areas, such as the United States and the 
Atlantic Ocean while the relation is poorest over the Pacific Ocean and 
North African areas. The 1.2 deg lat difference for Asia is considered 
anomalously low and not representative of that area. 

2. The objective HBZ positions north of 30N, and over eastern 
portions of continents and western portions of oceans, compare most 
favorably with the NMC positions. 

3. Only 66 percent of the NMC fronts could be associated with 
the objective fronts. The remainder are occlusions and/or weak fronts 
in eastern ocean areas. Both of these types are frequently analyzed 
from time continuity considerations only and are in the process of 
frontolyzing . 

4. A similar percentage of significant objectively-analyzed fronts 
(72 percent) could be associated with NMC fronts. The remainder represent 
baroclinic zones over North Africa and other subtropical ocean and land 
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areas as well as those fixed to orographic features. Following preva- 
lent air mass theory, the subjective analyst hesitates to place stationary 
fronts on the periphery of desert or mountainous areas, despite the nature 
of thermal patterns. 

A more recent, but quite limited, evaluation for the four synoptic 
times, 0000 GMT and 1200 GMT, 26 and 27 January 1967, considered approxi- 
mately 57,000 nautical miles of Northern Hemisphere fronts without regard 
to classification by region. An average position difference of 2.6 deg 
lat between NMC surface and objective (DGO^^^ model) 1000 mb fronts resulted. 
Forty-two percent of the numerically-computed fronts were found on the 
warm-air side of NMC fronts with an equal number on the cold-air side. 

The remainder were coincident, to within 1/4 deg lat. 

Tab_le 6. Statistical relations between numerically-analyzed 
1000 mb (DGG* model) and manually-analyzed NMC surface fronts for sixteen 
observation times in the period 23 August to 6 September 1965. 



Sector 


Objective Fronts Sampled 
(n mi) 


Average Distance 
Separating NMC & 
Objective Fronts 
(deg lat) 


North America 


49,250 


1.6 


Atlantic 


46,900 


1.4 


Pacific 


46,300 


2.1 


Europe-North Africa 


26,000 


1.8 


Asia 


18,000 


1.2 



Overall average 



1.7 
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850 mb Analyses 

Although, at this point, it is not the intent to discuss verti- 
cal consistency of numerically-analyzed fronts (see Chapter V ) it is 
appropriate to the substance of this chapter to show the 850 mb 
DG0* version for 0000 GMT 6 March 1966, in order to demonstrate the 
relative merits of numerical frontal analysis using mandatory leve 1 and 
layer data . 

From the conventional synoptic point of view, the 850 mb fronts 
are expected to be on the cold-air side of their 1000 mb counterparts. 

But, as noted from a comparison of Figs, 6 and 12, the objective 850 mb 
fronts are coincident with or are on the warm-air side of the 1000 mb 
fronts derived from the DG0* model. To show that this abnormal situ- 
ation is not uncommon, the following is offered in evidence. With 
reference to the evaluation of late January 1967 frontal analyses, 
described above, the average distance between 850 mb NMC and objectively- 
located fronts is 2.6 deg lat which is the same distance as that for the 
NMC surface versus objective 1000 mb. However, only 19,000 nautical 
miles of front were considered at 850 mb, most of which were over the 
dense -data North American region. Thus, considering the area sampled, 
the 850 mb positions of numerically-analyzed and NMC fronts compare less 
favorably than their counterparts at or near the surface. But, of greater 
importance is the fact that 86 percent of the 850 mb numerically-located 
positions were found on the warm side of NMC’s 850 mb fronts, with an alge- 
braic average distance of 2.3 deg lat. Such bias relative to NMC did not 
exist for the objective 1000 mb fronts derived from 0*. 
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Fig. 12. 850 mb DG0** analysis for 0000 GMT 6 March 1966. 

Isolines at interval of 20 x 10"2 c (100 km)^ , starting with +5 
units. The available NMC fronts at 850 mb are shown with conven- 
tional symbols. 
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Modification of DG0^ to Enhance the Analysis of Occlusions 

One of the more obvious omissions in the analyses generated by 
the DG0* model is the conventional occluded structure, as discussed 
earlier in the chapter* Occlusions are accepted and common to all 
national and international weather services, except perhaps that of the 
Canadians, whose working model features troughs of warm air aloft 
(trowals) as substitutes for the occlusion (Galloway, 1958). The well- 
developed occlusions on the surface charts of United States weather 
services are generally associated with a large-amplitude prefrontal 
tongue of warm air, as may be noted over the western At lantic--eas tern 
North America region and over the western North Pacific Ocean at 0000 
GMT 6 March 1966 (Fig. 7). Such occluded fronts intersect the isotherms 
at considerable angles and n^^ becomes a relatively poor approximation to 

—4 

n^ (See Eqs . 2 and 3). 

The close coincidence of the thermal axis of warm temperature and 
the conventional occlusion suggests application of the Laplacian operator 
to temperature in order to locate the **occluded** HBZ boundaries. Using 0* 
as the thermal field of interest, 

-V^G* = - V|V0*| *n-^ - |v0*|v-n-^ = DG0* + GK0*. (11) 

As noted in Eqs. (5) and (6), the DG0* parameter is indicative of the 
vorticity of the thermal wind due to shear and thus GK0* is representa- 
tive of the curvature component of the thermal vorticity. 

Fig. 6 demonstrates the poor relation of manually drawn occlusions 
to the field of 000*5^. GK0* (Fig, 13), which isolates the effect of the 
curvature term, is not adequate either to depict occlusions or any other 
frontal phenomena for that matter. The same comment is true far the 
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Fig. 13. GK0* (1000/700 mb) analysis for 0000 GMT 6 March 
1966. Isolines at intervals of 20 x 10"^ C (100 km)*^, starting 
with +5 units. NMC fronts are shown with conventional symbols. 
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Laplacian itself (not shown) since thermal tongues also occur in areas 
other than the location of conventional-front discontinuities. However, 
it is productive to combine the curvature and shear terms of Eq. (11) in a 
selective manner, that is only in the area of negative DG6* between the 
zero values within the HBZ and the adjacent subbaroc linic zone. This 
version is called Mod DG0* and is shown in Fig, 14, Such a procedure 
preserves the continuity of the non-occluded warm-side boundaries of the 
HBZs depicted by DG0* alone since these occur with positive DG0*. A blow- 
up of a section of the numerical frontal analysis over the western Atlantic- 
eastern United States area in Fig. 14 exemplifies the enhancement of the 
occlusion analysis by use of Mod DG0* (Fig, 15). Note that DG0* and Mod 
DG0* are the same except on the cold-air side of zero DG0* in the HBZ, 

It is of interest to compare Figs. 6 and 14. Note the differences 
in orientation and structure of the HBZ boundaries that are associated 
with occlusions, as for example over the extreme eastern and western 
sections of northern United States as well as northern Russia. Unfortun- 
ately, the curvature modification afforded by Mod DG0* is not always 
desirable. An example of this may be noted over the subtropical area 
of North America, a region frequently associated with low-tropospheric 
thermal tongues. Here, an already troublesome area for frontal analysis, 
due to topography and sparse-data problems, suffers further from the 
curvature modification. 

Overall, the major objection to the Laplacian adjustment for 
occluded structure is that the curvature-vorticity term is not directly 
related to baroclinicity and hence fictitious temperature gradients are 
implied. These objections have kept the program from becoming operationally 



useful. 
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Fig._14. (DG0* + GK0*) for DG0* < 0, GK0* > 0, and 

|gK0*| > 1dG0*|; 0000 GMT 6 March 1966. Isolines at intervals 

of 20 X 10"2 c (100 km)"^ starting with +5 units. Axis values > 
T5 units: xxxx : NMC fronts shown with conventional symbols. 
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Fig, 15. Modification to the DG0* (1000/700 mb) analysis for 
the case of a frontal occlusion over the western At lantic--eas tern 
United States area, 0000 6 March 1966. €* i^n K: light £olid 

lines; major axis of zero DG0*: dotted line; 000“* and mod DG0“* in 
10“2 C (100 km)"^ : heavy solid and dashed lines, respectively. 
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Some work has been done by Kirk (1965, 1966, and 1967) on incor- 
porating the Laplacian effect into objective frontal analyses. The 

2 

unqualified favorable appraisal of V T as a front- location parameter, 
given by Kirk in these references, completely misrepresents the utility 
of this approach. 

Effect of Mesh- length Variation on 
In Chapter II the effect of reduction in mesh length on the 
location and intensity of frontal zones was demonstrated for a limited 
case involving the modeled N and M temperature profiles. The results 
there suggest a similar study using real-atmosphere data. 

Consistent with both the simulation and operational experiments, 
the mean virtual potential temperature for the layer 1000/700 mb (O**) 
and the synoptic time, 0000 GMT 6 March 1966, were selected for this 
substudy relating mesh length variation to objective frontal analysis. 

In particular, the frontal zone chosen for study is that portion of the 
significant HBZ in the western At lant ic--eas tern North American region 
contained in the area identified by corners ABCD in Fig. 6. The subject 
HBZ is considered as moderate (M) in terms of the symbolic profiles in 
Chapter II (Fig. 3); the average baroc linic ity is 1.6C (100 km) ^ through 
the frontal zone. Starting with the full mesh grid-point values of 0*, 
the temperature data were numerically interpolated to the .75, .5, and 
.25 mesh intersections using Bessel*s central difference formula. Sub- 
sequently, the fields of and DG0* were produced for each of the new 

grids in a manner consistent with that described earlier in the Chapter. 

Fig. 16 illustrates the important aspects of the DG0* analysis 
on the full, half, and quarter mesh grids. Results for the three-quarter 
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Fig. 16. DG0* analysis as a function of mesh length for 
area identified by corners ABCD ^n Fig. 6, 0000 GMT 6 March 1966. 
Positive and negative axes of DG0*: from full mesh (heavy solid 

lines), half mesh (heavy^ long dashed lines), and quarter mesh (heavy 
short dashed lines). DG0* isolines intersecting positive (negative) 
axes, starting with +.35 (-.25) x 10“^ C (100 km)"^ at intervals of 
10 units: from full mesh (light solid lin^s), half mesh (light dashed 

lines), and quarter mesh (dotted lines). 0* isolines (K) indicated by 
edges of shaded areas. 
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mesh were intermediate to those for the full and half mesh and are not 
shown here. A background of full-grid 0* (shaded) is shown to point out 
an important departure from the modeled profiles in Chapter II, namely, 
that the warm air is much less barotropic and the cold air somewhat more 
baroclinic than the counterpart sections in the M profile (Fig. 3 and 
Table 1). However, at J < 15 similar conditions in the modeled and real- 
atmosphere temperatures do prevail. 

Analysis of the positions of warm- and cold-side boundaries of 
the HBZs and the frontal-zone intensities, as implied by the numerical 
DG9* analysis, indicates the following: 

1. Full versus half mesh: Overall, both full and half mesh 

patterns are similar in position and extent of the HBZ . Clear evidence 
of the positive correlation of mesh size and frontal-zone width is con- 
fined to the area below J14 where the latter is reduced up to 10 percent 
(40 kms) in the half- compared to full-mesh case. Poleward of J14 the 
half-mesh HBZ widths show both increases and decreases compared to full 
mesh, varying in amount from near zero to 25 percent. A large segment of 
the northern section shows both the warm- and cold-side HBZ boundaries 
displaced toward the warm air in the half-mesh compared to the full-mesh 
analys is . 

The reduction from full- to half-mesh relates well to the increase 
of intensity values, both along the warm- and cold-side boundaries of the 
HBZ, as expected from Chapter II. 

2. Full or half mesh versus quarter mesh: Fracturing of HBZ 

boundaries is very evident. In addition, the HBZ intensity does not 
uniformly increase from the .5 to .25 grid. Apparently, grid-size 
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reduction much below .5 is inconsistent with portraying synopt ic -sea le 
fronts from data originally fitted to the 380 km full-mesh grid. With 
an initial gross temperature resolution, the reduction to .25 of the 
full mesh (95 km) highlights mesoscale-like variations in the derivative 
fields. Such detail, partly derived from numerical-finite differencing 
and interpolation procedures, may be considered noise here since the data 
at full mesh, as a function of the mass structure model, do not contain 
such thermal resolution. 

It should be noted that some error in the relative positioning 
and magnitude of the numerically- analyzed fronts in Fig. 16 is due to 
the difficulty of locating the axes with an accuracy excelling .2 mesh 
length, even with the reasonably small analysis interval of 10 units of 
DG0*. 



Summary 

In summary, the results using real-atmosphere data are in general, 
but not specific, agreement with the simulated data experiment in Chapter 
II. It is evident that the resolution of the temperature data input to 
the DG0* program must be enhanced before reductions in grid-size much 
below 380 km can be effective in specifying the true width and intensity 
of frontal zones. Numerically- processed and analyzed thermal data, of 
the type utilized here, presently do not contain the detail inherent in 
data extracted from high quality subjective analyses. 

In order to focus attention on the problem of appropriate scale 
and resolution for acceptable numerical frontal analysis, the next 
chapter considers procedures for increasing the detail of the baroclinicity 
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data used in the DG0* computation. Highly resolved data, containing 
information on the mesoscale structure of the atmosphere, is readily 
extractable from vertical soundings of the atmosphere. 



CHAPTER IV 



EXPERIMENTS TO ENHANCE THE RESOLUTION OF BAROCLINICITY 



Most of the deficiencies of the operational DG0* model are due to 
the poor resolution of the temperature and baroc linic ity fields from which 
DG0* is computed. Frontal phenomenon, by its nature, tends toward a meso- 
scale structure while numerically-processed observations describe, at 
best, the synoptic -scale features. Therefore, some experimental effort 
was devoted to the problem of enhancing the resolution of HBZ analysis to 
a degree commensurate with the entity being analyzed. 

Modifications to DG0*, such as Mod DG0* in the last chapter, may 
improve the meteorological consistency of the present operational frontal 
model, but its merits cannot exceed the amount or scale of the information 
content inherent in the data tuned to a mesh size of 381 km. Thus, it is 
logical to look to the highly resolved significant- leve 1 radiosonde data 
for improving the HBZ analysis. One of the more promising experiments 
using detailed sounding data has its foundation in the analytical develop- 
ment which follows : 



0 = f (ng,p), 



( 12 ) 



where 0 is potential temperature; p is pressure; and n^, defined by 
^0 = V0/|V9| , represents an axis oriented in a direction perpendicular 



to an isentrope on an isobaric surface. Therefore, 

_ SG I ^6 r 



(13) 



6 
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Along an isentropic surface, d0 = 0, and from Eq . 13 

ro = ^ ^ 

5ng Sp 6ng ' 

Similar to Eq. (3), 

- 

DG0 = - VG0-n = - . 



( 14 ) 



(15) 



Computation of G0 and DG0 

The scheme proposed here involves a unique computation of G0 from 
the two factors, -S0/Sp and 6p/6nQ. This approach represents a consider- 
able deviation from that used to obtain G0* as a function of 0*, since G0 
is not derived directly from a field of 0. However, DG0 is computed from 
G0, similar to DG0* from G0*. 



The Computation 
of 50/5p 

The derivative, S0/Sp, is determined for a mandatory pressure sur- 
face by a quadratic fit of potential temperature data at the mandatory 
level and the first significant levels above and below. Plotted soundings 
on a thermodynamic chart give the impression of a sequence of layers, each 
linear in the logarithmic pressure ordinate, with separating first-order 
thermal discontinuities. In reality, the discontinuities bounding these 
significant layers are manufactured by criteria used for processing, 
coding, and plotting the observed data (Danielsen, 1959). Therefore, 
fitting a continuous curve to the transmitted sounding for the purpose 
of obtaining S0/Sp is reasonable as well as necessary. Also, using a 
quadratic rather than a linear fit of temperature data agrees with the 
nature of the sounding in pressure and temperature coordinates and ade- 
quately handles the situations in which the mandatory level is also a 



significant level. 
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The mandatory- leve L values of 50/Sp may be obtained through the 
solution for constants a, b, and c in the system of three simultaneous 
linear equations 



e 
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= a 


+ 


b 


Pu 
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= a 
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Pm 
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m 




= a 
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Pt 
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2 

cp^ , 



and from 



C|f) = + 2cp^, 



( 16 ) 



( 17 ) 



where. u(l) is the significant level geometrically above (below) the man- 
datory pressure level, m. The 0s in Eq . (16) are computed from sounding 
values of pressure and temperature by ^ 

I = 



In order to keep the program as simple as possible in the initial 
feasibility experiment for computing G0, temperatures to the nearest whole 
deg C (accurate to .5C) were employed in computing S0/Sp. This approach was 
partly dictated by the state of the hemispheric preprocessed radiosonde 
data made available by FNWC^. 

Significant layers of shallow depth (i,e., less than 10 mb) were 
not considered at this point in the experiment partly due to the gross 
temperature resolution and also because such layers may contribute intol- 
erable error in the S0/Sp computation. Further, as an example of the 
inconsequence of a thin layer, consider a stable significant layer near 



Since the data input to FNWC*s data-processing system excludes signifi- 
cant level information, it was necessary to write a complex program to 
couple significant and mandatory level data at the preprocessed stage. 
The program was written by Lieutenant R. G. Schiffner, Department of 
Meteorology and Oceanography, Naval Postgraduate School, Monterey, 
California . 
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700 mb of 10 mb thickness (approximately 400 feet) and 100 miles in 
horizonta 1 width . The slope of this layer, 6z/6n, is 1/1320. Hori- 
zontally wider and vertically thinner stable layers result in still 
smaller slopes and concomitant negligible baroc linicities along the 
pressure surface. Thus, HBZs of such dimensions may be profitably 
omitted from the frontal analysis, at least in this first attempt at 
computing G0, 

The Computation 
of 6p/5n 

Unlike SO/Bp, the computation of the slope of the isentropic 

surface, 6p/6n , at each station was accomplished within the framework 

of FNWC’s data processing and analysis system. Thus, 0* was substituted 

for 0 in Eq. (14). Available to the project were FNWC*s numerically-computed 

fields of pressure for each 0* surface at 2.5C intervals. From these 

fields, slopes of the isentropic surfaces, l(6p/6n) n |= 6p/6n , were com- 

P P 

puted by manual means. Here, n is the unit vector normal to isobars on 

P 

the isentropic surface. However , | (6 p/6n) = 6p/6n^^,the form of the slope 

factor in Eq. (14), represents only that component of 6p/6n^ which lies in 

the direction of V0* on the isobaric surface. Since the angles between the 
— > — * 

vectors n^ and n^^ were found to be < 15 deg for nearly all of the radio- 
sonde data used, the difference in the two quantities is negligible (i.e., 
generally < 3.5 percent in this particular experiment). In fact, over- 
estimating 6p/6n^ by using 6p/6n in Eq . (14) is considered desirable in 

0^ p 

view of the smoothing inherent in the numerically-processed products. 
Moreover^ 6p/6n^^ was computed for that isentropic surface whose magnitude 
is divisible by 2.5 and greater than 0* at the mandatory level for which the 
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computations are being made. For example, 6p/6n^^, for a station with 
0* = 301K at 850 mb, was computed at the 0* surface of 302. 5K. Undoubt- 
edly some fiction was introduced into the computed G0 by this simplified 
and hybrid combination of S0/Bp and 6p/6n^^. However, in the realm of a 
feasibility study such lack of homogeneity and detail may be tolerated, 
and, overall, may be desirable as will be discussed in a later section. 

Although the ultimate goal is to improve upon the operational 
DG0* model for frontal analyses, the emphasis here is on its predecessor, 
the baroclinicity , whose properties largely determine the validity and 
usefulness of the HBZ analysis. Thus, the discussion at hand will con- 
centrate on the enhancement of the field of baroclinicity as generated 
by the introduction of the vertical resolution inherent in significant- 
level radiosonde data. 

The gross scheme just outlined for computing G0 and its attendant 
field of DG0 will be symbolized as G0^ and DG0* to distinguish it from 
subsequent experiments which more closely approximate the true resolution 
of G0 and DG0. The test date is the same one used to exemplify the G0^ 
and DG0* computations, namely 0000 GMT 6 March 1966. 

G0* and its 
Comparison to G0^ 

In order to establish continuity with the discussion in Chapter 
III and present a base to which G0* may be compared, the numerically- 
analyzed 850 mb G0* for 0000 GMT 6 March 1966 (Fig. 17a) is discussed 
first. The analysis area is confined to North America and the middle 
and high latitude portions of the adjacent ocean areas from Europe to 
180 deg long. The positive DG0* axes closely parallel baroclinicity 
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magnitudes of 1C (100 km) , being on the weaker- (stronger-) gradient 
side for the less (more) intense HBZs . An example of the former (latter) 
may be noted in the middle (western) North Atlantic area. Baroclinic 
regions in which 9* gradients are less than 1C (100 km) ^ are of little 
interest but such gradients do occur along the DG0* axes. This is 
especially true for low-latitude ocean areas. Note the North Atlantic 
region south of 30N. The halfwidth of HBZ zones may be noted on Fig. 17a 
by the transverse distance from the DG0* to G0* axes. Average magnitude 
of the halfwidth is 2.5 deg lat. 

Fig. 17b shows the numerically computed G0' field, manually- 
analyzed in this version. About 40 percent of the North American area 
shows G0* less than 1C (100 km) ^ in the same area. Seven radiosonde 
stations yield G0* values greater than 3C (100 km) with an absolute 
magnitude of 8. 1C (100 km) ^ in southwest Canada, while only two G0* 
grid-point values exceed three units of baroclinicity . 

Table 7 statistically indicates the increased focus of baroclini- 
city using the significant-level radiosonde data from 118 stations in the 
area of analysis on Fig. 17b. The table shows the number of stations in 
each of nine categories of AG0 = G0* - G0*. The category **-.19 to +.19** 
may be considered as AG0 ~ 0. HBZs are arbitrarily defined as the areas 
having G0* values > 1C (100 km) For such areas, G0* > G0* in the vast 

majority of cases (85 percent), while in sub ba roclinic regions (SBZ) , 
where G0* < IC (100 km) G0* < G0* almost without exception (98 percent). 
This attests to the greater range of G0* relative to G0*, and especially 
to the reduction of baroclinicity in the s u b baroc linic areas on the G0* 
chart. Stated another way, the relative barotropic areas are more barotropic 
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while the relative baroclinic areas are more baroclinic in the G0’ as 
compared to the G0* fields. 



Table 7. 850 mb AG0 (= G0* - G0*) , 0000 GKT 6 March 1966, 

summarized for 118 radiosonde stations (North America and adjacent ocean 
areas), HBZ = stations where G0* > 1C (100 km)"^; SBZ = stations where 
G0*< 1C (100 km)“^. Stations with absolute instability omitted. 



AG© 


Number 


of radiosonde stations 


(C [100 km]-l) 


HBZ 


SBZ 


> + 2.00 


4 


- 


+1.00 to +1.99 


5 


- 


+.60 to +.99 


5 


1 


+.20 to +. 59 


14 


0 


-.19 to +.19 


13 


10 


-.20 to -.59 


4 


20 


-.60 to -.99 


1 


13 


-1.00 to -1.99 


0 


26 


< - 2.00 


0 


2 



The major axes have similar orientation on both the G0* and G0’ 
charts, but those on the latter are generally displaced to the cold side 
of their counterparts on the former. Note the situation over eastern 
United States and the adjacent Atlantic and the zone running northwest- 
southeast from Alaska to central United States, These relations between 
G0* and G0* (Fig, 12 and Table 7), suggest that the nature of the G0* 
computations reduces the broadening effect as compared to G0*, However, 
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an extra impetus to this effect may be due to the manner in which the 
slopes of isentropic surfaces (6p/6nQ^^) were computed. That is, using 
the ^higher 0* divisible by 2.5K’* for computing 6p/6nQ^ has the effect 
of moving the hyperbaroc linic zones toward the cold air from their true 
pos itions . 

It is also noted that detail along the G0* axes (eastern United 
States and southwest Canada) is greater than that of G0*. Secondary HBZs 
are both more evident and intense on the G0' analysis as seen on Fig. 17b 
over east central United States and southeast Canada. The mountainous 
regions are problem areas on both the G0* and G0* models; the two analyses 
differ markedly here. However, occluded structure, as shown by NMC over 
northwest United States and southwest Canada (Fig. 12), is indicated 
weakly with reasonable orientation in the G0* but not in the G0* fields. 

There still remain many problems in obtaining a suitable field of 
baroclinicity at any level, especially at low tropospheric levels. In the 
present experiment, the G0’ computation becomes invalid for the dry adia- 
batic case, that is for S0/Sp = 0 and 6p/6nQ^ = Also, the present 
mode of computing G0* is unrealistic for the occasional case of super- 
adiabatic conditions (S0/dp positive) in the observed sounding since this 
stability condition cannot exist in the numerically- processed data from 
which the bp/bn^. values were generated. Four such superadiabatic values 
occurred over the North American area for the synoptic time under experi- 
ment. This frequency would be greater, except that for such anomalous 
cases temperatures to the nearest .2C (vice 1C) were used in computing 
S0/5p. Teletype data transmissions were the source of such detailed 



temperature data. 
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DG0^ and its 
Comparison to DG0^^ 

To advance the experiment to the point of obtaining a DG0 * field 
for comparison to its counterpart, DG0*, the G0’ data were objectively 
analyzed. Several G0^ guess fields were considered. Even over a dense 
data area like the United States, a ”flat” guess, as G0* = 0, yields a 
numerical G0 ' analysis with a tendency for circular isolines due to 
localized points of large baroc linicity . Some improvement in shaping the 
field to elliptical rather than circular isolines is achieved by using 
the well-behaved G0'' field as first guess to the G0* analysis. Using 
this method. Fig. 18 shows a manual trace of the resulting 850 mb ob- 
jectively-analyzed G0' field. The reduction of resolution in G0 * -objective 
compared to G0 ’ -sub ject ive (Fig. 17b) is obvious. Also, differences in 
the positions of the axes of the two fields give a further measure of the 
role of objective analysis in misplacing the associated HBZ boundaries, 
a subject which came under discussion in Chapter II. 

The quantity DG0 ' (= -VG0*'n ) was computed next. The 9* field 

w ^ 

was used to generate the directional second derivative since a suitable 
numerical analysis of 0' was not readily available. 

Fig. 19 shows the DG0* and DG0 * numerical analyses, similar to 
Fig. 17 for G0* and G0*. In general, the comments concerning the baro- 
clinicity fields carry over to the interpretation of Fig. 19. Due to the 
close proximity of G0* baroclinic zones over the southeastern United 
States, the nearly meridional DG0* axis along the east coast is disrupted. 
This is an undesirable feature in view of the conventional-front positions 



for this time. 
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Also note that the axis of maximum baroc linicity in the G0* and 
G0* fields do not necessarily lie along their zero-isoline counterparts 
in the DG0’ and DG0^' patterns, respectively, as expected from the dis- 
cussion in Chapter III. The phenomenon is of significance only for com- 
plex HBZ structures which occur near the edges of dense data areas or in 
mountainous regions. The detailed G0*, DG0 ’ analyses are more susceptible 
than G0*, DG0* to this numerical problem, as noted over Mexico and the 
adjacent Gulf of Mexico (Fig. 19). 

Fig. 20 gives a measure of the positions of the HBZ boundaries 
from the DG0 ^ and DG0* models relative to subjectively-analyzed 850 mb 
NMC and Canadian fronts. DG0* excels in detail and in apparently reducing 
HBZ broadening by placing fronts on the cold-air side of the DG0* positions. 
However, the DG0^ fronts compare very favorably with the hand-analyzed 
versions in the matter of spatial continuity, especially over the western 
Atlantic Ocean. 

700 mb data and fields were handled in the same way as 850 mb and 
the results show similar merits and deficiencies. The vertical consist- 
ency of objective fronts from DG0* excels that from DG0*, as noted in 
Fig. 21. Improper frontal slopes (i.e., 700 mb positions on the warm-air 
side of 850 mb positions) are more prevalent in the DG0* than in the DG0’ 
patterns. The vertical consistency of DGO** fronts is illustrated and dis- 
cussed further in Chapter V. 

Before advancing the G0 experiment (Eqs . 12-18) in order to obtain 
a greater measure of accuracy and resolution in the baroclinicity fields, 
it is to be emphasized that the difference in detail between the G0 ' and 
G0* fields is derived from the computation of S0/Sp from significant-level 
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Fig. 20. 850 mb fronts, 0000 GMT 6 March 1966. NMC : conventional symbols; CAO, 

Canada; light solid lines (A = arctic, M = maritime, P = polar); DG9*: xxxx ; DG0' : oooo. 
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Fig. 21. HBZ boundaries (warm-air side 
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radiosonde date (Eq. 14). The influence of the 0^^ field is retained by 
substituting dp/6n^^ for 6p/6n^ in Eq. (14) and n^^^ for n^ in Eq. (15). 

Re computation of GQ with Certain Approximations Removed 
The previous section may be regarded as a successful feasibility 
experiment in obtaining a field of baroc linicity whose detailed structure 
excels that of due to the vertical resolution involved in the com- 

putation of S0/Sp. Utility of the derived G0* and DG0* is evident even 
though (1) error-prone unchecked real-time radiosonde transmissions were 
used to obtain S0/5p, (2) temperatures from such soundings were rounded 
to the nearest whole degree Centigrade, and (3) 6p/6n , at the mandatory 

level, was computed from numerically-processed data at a 0* surface whose 
value exceeded that at 850 mb by as much as 2.4C. 

Before proceeding to an extension of this experiment involving 
the enhancement of baroclinicity resolution by an increase in data 
density, it is well to remove several of the gross approximations in- 
herent in the computations to this point. For this purpose, the data 
for 850 mb, 0000 GMT 6 March 1966, were further manipulated. Due to the 
basic importance of G0 to the frontal parameter, DG0, only the former is 
in focus here. Concentration of effort is on the dense-data areas of the 
United States, Canada, and the surrounding ocean areas. 

For this second and improved approximation to G0 (hereafter 

symbolized as G0*) checked radiosonde data were used. The National 
c 

Weather Records Center, Asheville, North Carolina (NWRC, 1966) recomputes 
most soundings using both machine and manual processing. Among other 
things, the radiosonde data, as originally recorded , are checked for 



88 



proper evaluation and inconsistencies of mandatory- leve 1 data as evidenced 
by analyzed charts; also, derived height, temperature, moisture, and wind 
data are recomputed for 50 mb intervals as well as the significant levels. 
Such checking often results in soundings differing considerably from the 
transmitted version in both significant and mandatory temperature infor- 
mation. And, of course, both checked and transmitted soundings differ 
from the original recording and any systematic smoothing of same (Daniel- 
sen , 1959) . 

In view of the increased accuracy of the checked data compared to 
the transmitted soundings, the arbitrarily imposed accuracy of 1C in 
temperature and the 10 mb criterion for significant layers were removed. 
The following criteria were substituted in computing the baroclinicity 
parameter, . To compute S0/5p at 850 mb, temperature data were re- 
corded to the nearest .1C from 850 mb and the first significant levels 
above and below 850 mb. If such significant levels extended more than 
50 mb beyond the 850 mb level, then the 800 mb and 900 mb level data were 
used to formulate 50/^p. 

Use of the FNWC processed -data fields was maintained for computing 
the slope of the isentropic surfaces in Eq. (14). However, the slopes at 
each station were obtained from a more exacting differentiation scheme 
made available by Meteorology International (Roll, 1966). As noted in 
Chapters II and III, the derivatives used in the DG0* program were com- 
puted from the Lagrangian five-point differentiation formula for centered 
differences (Eq. [8]), which does not use the datum at the grid point of 

^The elaborate set of computerized checking procedures, called the 
Significant-Level Punched Card Audit, was made available to the writer 
by the NWRC. 
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computation. In Meteorology International's program, derivatives are 
computed from a 16-point 4x4 array of grid-point data with the com- 
putation point at its centrum. The program is a scheme which minimizes 
smoothing of the gradient values and weights heavily the influence of the 
immediate surrounding grid points as compared to the Lagrangian com- 
putations employed in the model. Further, for each radiosonde 

station, 6p/6n^. was interpolated for the 850 mb 9*, accurate to ,1C. 

A subjective analysis of the 850 mb data is shown in Fig, 

22. The pattern should be compared with those of G0' and G0* in Fig. 17. 

In comparison with G0* , G0^ shows only a slight (about 10 percent) reduction 
in the subbaroclinic area but the stations having baroclinicity magnitudes 
> 3C (100 km) ^ have increased to twelve in number. These figures agree 
in principle with a sampling of data from those radiosonde stations for 
which there exists both G0' and G0^ values; namely, 63 percent of the HBZ 
baroclinicities have increased while only 36 percent of the SBZ values 
have decreased from G0* to G0'. The absolute high magnitude of G0* is 

C o o ^ 

21,30 C (100 km) ^ at station 879 in southwest Canada, thus increasing 
considerably the range of G0^ relative to G0*. 

Similar to G0* versus G0* (Table 7), G0^ gives more focus or 
resolution than G0* (Table 8) with 82 percent of the HBZ values increasing 
and 89 percent of the SBZ values decreasing from the G0* to the G0^ com- 
putations, It is to be noted that there is a 25 percent increase in 
comparative data points in Table 8 over that in Table 7, 

There are some important differences in the patterns of G0' and 
G0^. The most significant of these occur over the southwest and southeast 
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United States, the former associated with orographic effects, the latter 
with a synoptic-scale frontal structure whose complexity is exemplified 
by the Canadian analysis in Fig. 20. Repositioning of the axis of maxi- 
mum baroclinicity in eastern and central Canada is also apparent. There 
can be little question from Fig. 22 that the accuracy and resolution of 
the radiosonde data are important factors in locating the axes of baro- 
clinicity and hence the hyperbaroc linic zones. 

Table 8. 850 mb AG0*(= GQ^ - G0*), 0000 GMT 6 March 1966, sum- 

marized for 147 radiosonde stations (North America and adjacent ocean 
areas). HBZ = stations where G0^ > 1C (100 km)’^; SBZ = stations where 
G0^ < 1C (100 km)’^. Stations with absolute instability omitted. 



AG0 * Number of radiosonde stations 



(C [100 km]"^) 


HBZ 


SBZ 


> + 2.00 


9 




+1 . 00 to +1 . 99 


6 


- 


+.60 to +.99 


9 


0 


+.20 to +. 59 


17 


5 


-.19 to +.19 


16 


34 


-.20 to -.59 


7 


23 


- . 60 to - . 99 


1 


6 


-1.00 to -1.99 


1 


12 


< -2.00 


0 


1 



In general, G0^ differs from G0* by virtue of greater accuracy and 
detail in computed BO/Sp and a reduction of smoothing and increase in 
accuracy of 
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Many problems still remain. Nine radiosonde stations, labeled 
on Fig. 22, indicate absolute instability in the layer within 50 mb 
of 850 mb. Thus 6p/6n^^^, as derived from FNWC’s mass structure model, 
is in error at these stations. A potentially greater problem, which 
encompasses the absolute instability difficulty, is the validity of 
6p/6n^^ in combination with SO/Sp. 

Cross-section 
Ana lyses 

A cross-section analysis, extending from station 78367 (Guan- 
tanamo Bay, Cuba) to 72655 (St. Cloud, Minnesota) is used as another 
means of exemplifying the relation of baroclinicity from numerically- 
processed and checked radiosonde data. Fig. 23 shows the former, with 
0* data plotted at 50 mb intervals, while Fig. 24 shows the latter, with 
0 computed from NWRC’s Northern Hemisphere data tabulations (1966) using 
Eq. (18). In the case of station 78367, significant- leve 1 data were 
taken from the real-time transmitted sounding, as checked data were not 
available. The mixture of checked and unchecked data is justified for 
this station since mandatory level data on the two versions are nearly 
identical, a fact not true for all radiosonde stations. 

The subject cross-sections, oriented approximately in the 
direction of V0* at 850 mb, are quite interesting for several reasons. 
The sections intersect a complex frontal structure over the Florida 
area (see Figs. 20 and 22), extend from a near maximum to a near minimum 
in 850 mb temperatures (Fig. 10), and pass through a relatively dense 
data area . 

Contrary to the expectation that 0* should exceed 0, nearly 
40 percent of the area from the surface to 600 mb shows 0 > 0*, in some 
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Fig. 23, Vertical cross-section from 78367 to 72655 at 0000 GMT 6 March 1966. 0* (K) 

plotted to left, computed from reconstructed soundings of FNWC numerically-processed mandatory 
level data. Isoline interval: 5 K. 0* - 0 (K) plotted to right; (0* - 0) indicates values 

estimated from analysis. 
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Fig. 24. Vertical cross-section from 78367 to 72655, 0000 GMT 6 March 1966. 0 ( 
plotted to left, and T (C) plotted to right, computed from checked radiosonde data (except 
72367). 0 isoline interval: 5 K. Stable layer (-dT/Sz < 0) boundaries: dashed lines. 



96 



cases (station 74794, at 1000 mb) by as much as 5C, The situation is 
similar at the 850 mb level. Generally, 0* exceeds 0 in the relatively 
warm, moist areas adjacent to the frontal zones, while in other areas 0 
is greater than 0*. To assist in this analysis the quantity 0* ~ 0 is 
plotted at 50 mb intervals in Fig. 23. 

Although the major axis of warm and cold air have similar locations 

on the two versions, detail on stability (from station 78367 at lower levels 

to 72327 at 650 mb) and positions of the minor thermal axes differ markedly 

in pattern (as 900-650 mb, from 78367 to 72206). It appears that isen- 

tropic slopes on the numerical version are excessive in stable areas 

(which would lead to excessive G0* values in HBZs) and deficient in or 

c 

near actual unstable areas (leading to deficient G0^ values in subbaro- 
clinic and unstable layers), relative to slopes derived from the checked 
radiosonde data. The former is to be noted in the vicinity of the sloping 
stable layer (s) equatorward of station 72311 while the latter is obvious 
at the near-surface levels throughout the cross section. The situation 
at low levels (notably at 72327 and 72655) also serves as an example of 
the inability of the numerica 1-data-processing scheme to specify the 
location of absolutely unstable layers in the actual atmosphere. 

Broadening effects due to the numerical processing of data are 
quite obvious in Fig. 23 within and either side of the major zones of 
hyperstability and hyperbaroclinicity , especially from 850 mb to 700 mb 
between 78063 and 72311. However, the reverse is somewhat true geo- 
metrically below 900 mb from just south of 78063 to 74794, where 0* 
exceeds 0 (moisture effect) while 0* is generally less than 0 north of 
78063. This is a propitious circumstance for locating the HBZ in the 



numerical field of data. 
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It is of interest to seek further quantitative comparison of the 
isentropic slope values since the mass -s tructure version of this factor, 
6p/6ng^, is retained in computing Table 9 shows the calculated 

850mb slopesateach station as derived from the two cross-sections as well 
as those generated from the numerically-processed data fields 

and used to compute G0^ in Fig. 22. 

Measurements of 6p/6n and Baroclinicity --85Q mb point slopes of 
isentropic surfaces from the cross sections of analyzed processed data 
(Fig. 23) and checked radiosonde data (Fig. 24), and average slopes used 
in computing G0* (Fig. 17b) and G0^ (Fig. 22), are compared for the six 
stations from 78063 to 72532 in Table 9. The range of isentropic slopes 
is greatest from the point slope computations made from the cross- 
sections. Beyond that, clear cut differences do not exist between the 
various values of 6p/6n. 

The 6p/6n values may be combined with 30/5p computed from the 
checked data to produce the baroclinicit ies , G0, shown in Table 9. The 
pattern of G0s associated with the isentropic point slope measurements 
from the cross-section in Fig. 24 agree with the complex doublet of HBZs 
between stations 72206 and 78367. Such is not indicated in the other 
versions of G0 in Table 9. Somewhat disturbing is the relative maximum 
of G0 in the vicinity of 72327 for all four G0s . Since this feature is 
largely due to the steep isentropic slope, the G0 version which reduces 
its significance is most desirable. Overall, as expected, real data vice 
processed data and point-slope vice average-slope calculations yield a 
G0 picture closest to the mesoscale nature of HBZs. 
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Table 9. 850 mb 6p/6n (mbflOO km]’^)0000 GMT 6 March 1966: 

Computations of point and average values for stations 78063 to 72532 
in cross-sections of Figs. 23 and 24. G0*s(C[lOO km] ; computed from 
tabular values of 6p/6n and SO/Sp computed from checked radiosonde data. 



Station 


Point slope 
from cross-section 


Average slope numerically 
computed from 




Checked data 
(Fig. 24) 


Processed data Met. Inti. 

(Fig. 23) Program 


Lagrangian 
5-point 
for mu la 








6p/ 6n 






78063 


48.5 




28.5 


41.0 


28.0 


74794 


21.0 




44.0 


49.0 


-- 


72206 


27.0 




62.0 


45.0 


35.0 


72311 


104.0 




90.0 


75.0 


52.0 


72327 


64 . 0 




150.0 


60.0 


72.0 


72532 


36.0 




25.5 


36.0 


20.0 








G0 






78063 


3.1 




1.8 


2.6 


1.8 


74794 


1.0 




2.2 


2.4 


-- 


72206 


1.4 




3.2 


2.3 


1.8 


72311 


0.9 




0.8 


0.6 


0.4 


72327 


1.0 




2.2 


0.9 


1.1 


72532 


0.4 




0.3 


0.4 


0.2 


The 


resolution 


of baroclinicity is 


advanced one more step in the 


next section 


where the 


detail 


and accuracy 


of the analysis 


is increased 



by a scheme of extrapolating significant-layer G0 data to the 850 mb level 
at points between the locations of radiosonde stations. 
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G0 Extrapolations 

The increased resolution of the field of as compared to G0* , 

is due to the computation scheme and nature of the data; the data density 
is essentially the same for both. The average distance between radiosonde 
stations is close to 5 deg lat over dense data areas and much larger than 
this over ocean areas, thus limiting the enhancement of detail in the 
analysis of baroc linicity for any G0 computation scheme. A partial 
solution to the data problem is made possible by an extrapolation technique 
which has its basis in the thermodynamics and synoptics of the vertical 
structure of HBZs . It is well recognized that a hyperbaroclinic zone of 
significance is generally associated with a relatively stable layer of 
air in which S0/Sp is negative and large. Vertical cross-sections, 
oriented approximately along V0 on an isobaric surface, show such layers 
extending over atmospheric depths of several thousand feet and distances 
of hundreds of miles with boundaries identified by an isentrope or a small 
range of isentropes (Saucier, 1955). That stability and, consequently, 
baroclinicity have considerable uniformity in a particular stable layer 
may be seen from the cross section previously discussed (Fig. 24). 

Based on this general background, selected layers above and below 
a mandatory pressure surface have been examined for stability (50/Sp) and 
isentropic slope (6p/6n) for the purpose of extrapolating baroclinicity 
data to the mandatory level at locations between the radiosonde stations. 

In this way, the vertically-resolved temperature data contribute further 
to the detail of the baroclinicity analysis. The specifics of the extrapo- 
lation procedure finally adopted are outlined below. Application is made 
to the 850 mb level, 0000 GMT 6 March 1966. 
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Sq/5p - -E xtrapolations are derived from radiosonde data between 
700 and 950 mb only. So/Sp is computed at 900, 800, and 750 mb, using 
a quadratic fit of checked radiosonde temperature data. For each of the 
three levels, data points are selected for the level and 50 mb below and 
above unless an intervening significant level is available. For example, 
to compute So/^p at 900 mb, 0s at 950, 900, and 850 mb are used if signi- 
ficant levels are absent between 950 and 850 mb. If significant levels 
existed at 920 and 855 mb, then data at the 920, 900, and 855 mb levels 
are used in the So/Sp computation. Such a procedure obtains three 
extrapolations for each radiosonde station (unless station pressure 
< 950 mb) and controls the layer from which the extrapolations are 
derived (i.e., from approximately 2500 feet below to 5000 feet above 
850 mb). Due to the complexity of localized stable layers near the 
ground, the use of data in the low troposphere below the 850 mb level 
must be limited. A less arbitrary system could set the lower-level 
limit at 50 mb above the earth’s surface, vice a fixed 950 mb, thereby 
handling data in mountain and flat areas in similar fashion. 

5p/ 6 n - -Similar to the previous application of Eq. (14), this param- 
eter is computed as 5p/6n^^ at 900, 800, and 750 mb. 

G9”--The product of S0/Sp and 6p/6n^^ is S 3 nnbolized as G0” for 
the cases of data computed at P ^ 850 mb. The computed G0” values are 
extrapolated to 850 mb along the 0* surface intersecting the point of 
computation at pressure level P (i.e., 900, 800, and 750 mb) and in a 
direction along the 850 mb V0*. 

It is to be noted that the system adopted here extrapolates 
baroclinicity to 850 mb rather than simply extrapolating stability. 
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So/Sp, along the isentropic surface. One reason is that the accuracy 
of , using a 6p/6n^^^ computed at P ^ 850 mb for the known position 
of the radiosonde station, is better than that obtained from use of 
6p/6n at an estimated position of 0* at 850 mb, 

A considerable amount of manual manipulation of data was necessary 
in the extrapolation phase of the project since this aspect was handled 
as a feasibility study. Although S0/Sp was numerically computed at 
P ^ 850 mb from checked radiosonde data, the 0* appropriate to the 
6p/6n^. computation was determined manually from P, 0* soundings 
reconstructed from numerical grid prints of P (0*) at 2.5K intervals. 
Accuracy of interpolated 0* is considered to be + ,3K, For values of 
0* at 900, 800, and 750 mb^the isentropic slopes were interpolated. 

Analysis of Combined 

G0* and G0" Data 
c c 

Fig. 25 represents an analysis of the combined baroclinicity 
extrapolations from 900, 800, and 750 mb to 850 mb and the G0^ compu- 
tations at 850 mb, the latter analyzed in Fig. 22. The extrapolations 
are plotted to the nearest .1 C (100 km) ^ to distinguish them from the 
850 mb data plotted to the nearest .01 C (100 km) Since the atmos- 
phere is constrained to be statically stable in the numerically-processed 
data (Chapter III) used for 6p/6n^^,^, the 900 (800 and 750) mb extrapo- 
lations always appear on the cold-air (warm-air) side of each radiosonde 
station. For a few stations, as 72365, 72485, and 72576, extrapolation 
data only exist since station pressures are less than 850 mb. The 
extrapolation data are particularly helpful in areas where the algebraic 
sign of 850 mb S0/Sp is positive (absolute instability cases, plotted as U) , 



as for stations 72232 and 72374, 
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Fig. 25. 850 mb G0’ and G0” extrapolations, 0000 GMT 6 March 

1966. Analysis interval: geometric, with r = 2; units: C (100 km)"l. 
Axes: dashed lines. 850 mb G0^ data (radiosonde-station circles 
shaded) plotted to nearest .01; extrapolation data, G0”, (location 
circles not shaded) plotted to nearest .1. U indicates S0/Sp is 
positive; 850 mb and extrapolated data for each station connected 
by light dashed line. 
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Slopes of the isentropic surfaces (6p/6ng^) are indicated by the 
distances of the G0^ extrapolations from the parent radiosonde station. 
These distances correlate inversely with the magnitude of As an 

example, compare extrapolations for 78016 and 72206; the latter (former) 
indicates nearly vertical (horizontal) isentropes and relatively high 
(low) baroclinicities , Validity of a baroclinicity computation decreases 
with increasing distance from the radiosonde station (i.e., decreasing 
isentropic s lope) . 

Extrapolations were not plotted for stations where the 0* surface 
at P ^ 850 mb could not be located at the 850 mb level or where V0* at 
850 mb changed by 180 degrees between the station and its associated 
extrapolation point. For example, this situation occurred at 72867 
(76255) for the 900 (750) mb extrapolation. A 0* of 268. OK (309.8K) 
did not exist at 850 mb. For this study all the extrapolations (except 
those indicated with a question mark) were held to be valid and used in 
the analysis. The resulting mesoscale detail in Fig. 25 may be questioned. 
Perhaps those extrapolation points separated from the parent radiosonde 
station by a baroclinicity axis should be given minimal weight in the 
analysis. The same approach may be appropriate for data associated with 
nearly horizontal isentropic surfaces and isolated points of excessive 
high and low G0” values. Subjective computations and use of the hybrid 
mixture of 0* and 0 information very likely inflates the baroclinicity 
magnitudes in the relatively stable layers but deflates these values in 
the layers with instability, as noted from a discussion of Figs, 23 and 



24 . 
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In general, the number of baroc linicity axes from G0^* and G0^ 
data combined have doubled compared to G0^ alone. This relation between 
the two fields is much like that occurring in a comparison of flight 
and radiosonde data for the purpose of locating the jet-stream axis. 

Each synoptic-scale jet axis may be reconstructed into a multiple number 
of splinter axes by increased resolution of the observational data. 

Although the analyses over the eastern United States-wes tern 
North Atlantic region appears complicated it is similar in interpretation 
to the detailed Canadian CAO frontal analysis shown in Fig. 20. For 
instance, the CAO analysis shows 850 mb arctic, maritime, and polar fronts 
distinguishable over the Florida Peninsula. These frontal entities are 
in much the same position on Fig. 25. For the remainder of eastern 
North America, the comparison in number of CAO and numerical HBZs is 
good, but the positions of these zones differ. 

The profiles of 850 mb baroc linicity between stations 78367 and 
72655 on Figs. 24 and 25 are quite comparable, except near 78367. A 
complex system of hyperbaroc linicity between stations 78063 and 72206 is 
found on both charts as is the axis of maximal values near 72655. Station 
72532 presents a problem in correctly extrapolating data to the 850 mb 
level since an isentropic minimum occurs in the layer 800-700 mb very 
close to this station. 

Over mountain areas the analysis from G0’ and G0” combined differs 
sharply from that of G0^ or G0* alone. In fact, differences between any 
two of the manually or numerically generated baroclinicity analyses are 
maximized over such areas. A satisfactory solution to meaningful patterns 
over these regions has yet to be achieved. 



105 



Summary 

The various experiments to enhance the resolution of baro- 
clinicity for the purpose of HBZ analysis suggest the following: 

1. There is considerable baroc linicity detail inherent in the 
vertically-resolved significant- leve 1 radiosonde data. Such mesoscale 
detail cannot be retrieved from numerically-processed mandatory- leve 1 
data alone. 

2. Since significant- leve 1 radiosonde data (transmitted or 
checked) are of a meso-microscale nature, the baroc linicity values 
derived from them may not be compatible with the horizontally-projected 
synoptic-scale fields of other meteorological parameters. However, 
numerically tuning the baroc linicity field by use of smoothers and 
filters can solve this problem to any desirable degree. 

3. The numerical methods and grids in current operational use 
for processing and analyzing meteorological data are not appropriate to 
the scale of baroclinicity shown in Fig. 25. For this reason, the develop- 
ment described in this chapter must be regarded as a second-generation 
model for the purpose of analyzing hyperbaroclinicity . 

Next, attention is returned to the mode of HBZ analysis feasible 
for present-day field operations, which is described as the DG0* model 
in Chapter III. Considerable insight on synoptic-scale structure and 
behavior of fronts may be gained from a diagnosis of this model in 
varying space and time frames and in relation to other meteorological 
parameters. Thus, Chapter V concentrates on the four -dime ns iona 1 con- 
tinuity of HBZs , as illustrated by case studies from winter and summer 
periods. Chapter VI relates the HBZ analyses to weather satellite 
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observations of cloud and vertical motion; and, Chapter VII demonstrates 
the properties of time-averaged intensity and location of HBZs for a 
selected 30-day period in both winter and summer. 

Together, the next three chapters give the synoptic meteorologist 
a four-dimensional view of frontal analysis not heretofore obtainable 
from the standard manual analysis of surface fronts only. 



CHAPTER V 



THE DG0^ MODEL: CASE STUDIES OF THE SPACE AND TIME 

CONTINUITY OF HYPERBAROCLINIC ZONES 

A Summer “Sea son Case Study 

The summer-season day of 20 August 1965, with emphasis on the 
0000 GMT synoptic time, was selected primarily for a study of the space 
continuity of the tropospheric DG0* fields. Analyses at 1000, 850, 700, 
and 500 mb were considered in the substudy, 

DG0* Analyses at 1000, 

850, 700, and 500 mb 

The DG0* operational model developed in this treatise explicitly 
indicates boundaries and implicitly the intensities of HBZs , In general, 
most meteorologists are familiar with the synoptic-scale boundaries and 
their vertical-plane slopes from an air-mass point of view. For instance, 
it is well known that, in the low-troposphere , frontal zones slope toward 
cold air with increasing height, with slopes generally steepest at the 
lowest levels. In the middle troposphere, these zones of hyperbaroclinicity 
become more horizontal and diffuse. However, not much is known of the 
intensity of HBZs, since, normally, the parameter is not displayed in a 
quantitative fashion on the daily operational charts. 

It is the purpose of this case study to diagnose the numerically- 
computed boundaries of HBZs as generated from the processed virtual 
potential temperatures (0*) and to relate these analyses to the previous 
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discussion concerning the effect of numerical processing on the location 
and intensity of HBZs (Chapter II). 

Figs. 26-30 illustrate the DG0* fields and related phenomena for 
0000 GMT. 20 August 1965. Isolines on these figures represent a manual 
trace of selected positive- and negative-va lued lines of the numerically- 
produced front-location parameter. 

In order to focus attention on the HBZs with spatial significance 
and continuity, certain criteria involving intensity and pattern were 
established before locating the axes in the DG0* fields. The major 
criteria are as follows. (1) Analysis is limited to the area poleward 

j ! 

of 15N. (2) Closed isolines are not drawn around single grid-point values, 
whatever the magnitude of DG0*. (3) The positive and negative axes are 
drawn only for segments along which both associated axes are well defined; 
DG0* must exceed 10 units at one or more grid points along one of the axes. 

^ ■ M r 

(4) Axes are drawn only where DG0* magnitudes exceed 5 units. However, 
such axes are connected through a DG0* area of less than 5 units magnitude 
(but greater than zero) where the interval of separation is less than 2 
mesh lengths, and the associated positive or negative axis is continuous 
using the 5-units criterion. See DG0* analysis over northern United 
States on Fig. 26. (5) Axes associated with areas of circular DG0* isolines 

are not drawn. See Kamchatka, northern Russia, and the southeast China 
coastal areas in Fig. 26. Such areas generally signify cut-off pools of 
cold air or complex isotherm patterns. (6) Unless a HBZ is associated 
with at least one grid-point G0* value in excess of 1 C (100 km) ^ the 
DG0* positive and negative axes are not shown. Such is the case for many 
ocean and low-latitude HBZs. (7) Furthermore, axes obeying the criteria 
in (1) to (6) are also required to be approximately parallel although the 
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transverse distance separating the HBZ boundaries is of no particular 
limitation . 

Establishing criteria of this character has as its purpose an 
objective separation of the significant from the insignificant synoptic- 
scale HBZs . This selective screening may be also accomplished on an 
operational basis. Frontal zones weaker than those in focus here are 
considered as noise; the complex zones omitted by application of 
criteria (1) to (7) generally involve undesirable fiction due to the 
numerical processing in combination with grid size and/or scale of the 
thermal feature. 

1000 mb (Fig. 26) --The HBZs at 1000 mbs show less horizontal 
continuity and have more complex structure than their free -atmos phere 
counterparts. Features existing at 1000 mb but not found at success- 
ively higher levels, as over mountainous terrain, are of little real 
consequence. In general, oceanic positive and negative axes are parti- 
cularly sporadic, especially so over the Pacific Ocean. The DG0* 
intensities are greatest over land. The pronounced HBZ over western 
North America, a persistent summer feature, is a function of the 
differential large scale heating either side of the land-sea interface, 
with some contribution from the regional trade wind inversion (Neiburger 
et al., 1961). Occlusions over land and ocean areas are poorly indicated 
as expected from discussions in Chapters II and III. 

850 mb (Fig. 27 ) --Smoothing and reduction of segmented HBZs are 
evident when compared to 1000 mb. Intensity of baroc linicity change at 
the cold boundaries of the HBZs appears less than at warm boundaries for 



the non-s tationary frontal zones over east continents -wes t oceans, as 
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Fig. 26. DG0*, 1000 mb, 

10 X 10"2 c (100 km)"2, starting 
isolines: solid (dashed) lines. 

HBZs : heavy solid (dashed) lines 



0000 GMT 20 August 1965. Interval: 
with the 5 isoline. + (-) DG0* 
Warm-air (cold-air) boundaries of 
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Fig. 27,- DG0*, 850 mb, 0000 GMT 2C* August 1965. Interval: 

10 X 10**2 c (100 km)"2, starting with the 5 isoline. + (-) DG0* 
isolines: solid (dashed) line;s . Warm-air (cold-air) boundaries of 

HBZs : heavy solid (dashed) li]ies. 




112 



expected from the analysis in Chapter II. The quas i-station iry pattern 
due to differential heating over western North America is st.ll evident 
at 850 mb . 

7QQ mb (Fig. 28) -"More extensive and continuous fronzs 'prevail 
at 700 mb compared to lower levels. Reorientation of baroc 1 Lnicity on 
the North American west coast indicates the surface- induced thermal 
structure does not exist at this level. HBZ boundary values generally 
show an increase from lower levels, especially over ocean areas. 

500 mb (Fig. 29) --Operational use of frontal zones at levels 
geometrically above 700 mb is practically non-existent due to the con- 
centration on low-level air-mass structure. However, 500 mb frontal 
phenomena do exist, as evident from Fig. 29. This level shows 
latitudinally extensive and intense HBZs which take on the appearance 
of spiraling hemispheric-scale jet axes (Riehl, Alaka, Jordan, and 
Renard, 1954). Many of these zones do not have apparent counterparts at 
the lower levels. Previous authors have noted a similar relation between 
the synoptic-scale HBZs of the lower and upper troposphere, as well as the 
relative importance of the latter to movement and development of tropos- 
pheric weather systems (Sanders, 1967). 

Vertical Continuity 
of HBZ Boundaries 

The positions of the warm-side boundaries of the HBZs at the five 
levels, 1000, 850, 700, 500, and 300 mb, have been included here to empha- 
size the status of vertical continuity in the operational DG0* model (Fig. 30). 
The frontal contours that lie over the eastern United States, Canada, the 
northern Atlantic, and northwestern Europe areas are internally consistent. 
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Fig. 28. DG0* 700 mb, 0000 GMT 20 August 1965. Interval: 
10 X 10"^ C (100 starting with the 5 isoline. + (-) DG0* 

isolines: solid (dashed) lines. Warm-air (cold-air) boundaries of 

HBZs : heavy solid (dashed) lines. 
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Fig. 29. DG©*, 500 mb, 0000 GMT 20 August 1965. Interval; 

10 X 10"2 C (100 km)“2, starting with the 5 isoline. + (-) DG0* 
isolines: solid (dashed) lines. Warm-air (cold-air) boundaries of 

wR7c • h#:»;i'V7v qnlid Tdashed') lines. 
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Fig, 30. 
20 August 1965. 
/// : 300 mb : xxx 



Vertical continuity of positive DG0* axes, 0000 GMT 
1000 mb: ; 850 mb: ; 700 mb: **•; 500 mb: 



116 



Thus 5 slopes appear most logical for those frontal zones which are not 
tied to orographic features and where data are ample. Removal of 1000 mb 
fronts would result in further improvement in vertical consistency, 
particularly over Asia. The weak low-level fronts over oceans have 
poorest spatial continuity, as manifested by HBZs tilting toward warm 
air with increasing height. With reference to Chapter II, improper 
slopes may be expected for HBZs which intensify with increasing height, 
since the relative broadening effect increases with the frontal-zone 
strength. However, distinctly separate upper-tropos phere fronts may 
account for the seemingly poor continuity between upper and lower 
troposphere over the ocean areas. Moreover, the liklihood of tropo- 
pause intersections with the 300 mb surface further complicates the 
continuity shown in Fig. 30. 

Although not shown here, the front- locat ion parameter was 
computed from the objectively-analyzed field of unprocessed observed 
dry-bulb temperatures at 1000, 850, 700, and 500 mb, 0000 GMT 20 August 
1965. Maximum differences existed between DG0* and the dry-bulb version 
at 1000 mb. Negligible differences existed at 500 mb and 700 mb. It 
was also evident that the smoothing inherent in numerically-processed 0* 
somewhat compensates for the potentially greater range of 0^ compared 
to 0 at a 11 leve Is . 



A Winter-Season Case Study 

This section considers vertical consistency and time continuity 
of numerical frontal analysis for the winter period 1200 GMT 25 January 
to 0000 GMT 28 January 1967. As before, 1000, 850, 700, and 500 mb 
levels are subjected to the diagnosis. In this case the area is limited 
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to the North American section of the Northern Hemisphere equatorward 
of 55N latitude, in order that conclusions be drawn from evidence in 
dense-data areas only. 

HBZ Boundaries in 
Vertical Cross Section 

The medium of presentation used here is largely that of the 
vertical cross section depicting HBZ boundaries in the layer 1000 to 
500 mb, as derived from mandatory pressure- leve 1 analyses of 0* and 
DG0* information. For this series, consideration was given only to 

-2 

those HBZs marked by six or more grid-point values of DG0* > 5 x 10 

_2 

C (100 km)“ , at least one of which > 15 units. 

Fig. 31 shows a schematic 1000 mb frontal contour as an aid in 
identifying the locations of the cross-section groups in Figs. 32-34, 
Nineteen vertical cross sections, intersecting areas adjacent to develop- 
ing frontal-wave cyclones over central and eastern United States, were 
selected for analysis. The sections are categorized into groups; each 
group is homogeneous relative to the position of the apex of a frontal- 
wave cyclone and the movement of the surface front. Both cold- and 
warm-side boundaries of the HBZs are shown in each case. Boundaries 
tilting to the left (i.e., toward cold air) with increasing height have 
positive or normal slope, 0* values at HBZ boundaries, for each mandatory 
level considered, are also shown. Some of the quantitative conclusions 
which may be drawn from this sample of nineteen cross sections are shown 
in Table 10 and discussed in the following paragraphs. 

HBZ Width - -Sawyer (1955), interpreting British flight data taken 
at levels ranging from 1,600 to 20,000 feet, indicates that the widths of 
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frontal zones average 115 n mi, with a contrast of 5C(i.e., 2,5 C 
[deg lat] ^), closely matching the N profile in Fig. 3, However, if 
radiosonde data are used to determine the air-mass contrasts for the 
same times as the twenty observation flights then the average width 
becomes 525 n mi, with a contrast of 8 C (i,e.,,9 C [deg lat] ^) , 
similar to the W profile. Expressed as a range. Sawyer’s frontal-zone 
widths vary from 50 to 225 n mi. More recently, Danielsen has shown 
analyses of frontal zones over the western United States, derived partly 
from middle- and high-troposphere flight data in support of Project 
Springfield (Danielsen, 1964). Baroclinicity magnitudes computed from 
Danielsen ’s selectively-sampled stable zones above 2,000 feet fall in 
the range 12 C (80 n mi) ^ (i.e., 9 C [deg lat] ^) to 1.5 C (40 n mi) ^ 
(i.e., 2.2 C [deg lat] ^) . Such flight-data values generally exceed those 
of the hypothetical profiles of Fig. 3, which were designed to be typical 
of baroc linicities diagnosed from a constant-pressure analysis of radio- 
sonde data^as indicated below. 

The numerically-analyzed HBZ widths in the January 1967 cross 
sections average close to 400 n mi, increasing slightly with increasing 
height. Average HBZ intensity decreases with increasing height from 
3.5 C (deg lat) ^ at 1000 mb to 1.6 C (deg lat) ^ at 500 mb, thus ranging 
from the N to M profiles in Fig. 3. As expected, most of the associated 
DG0* values at both the warm-side and cold-side boundaries of the HBZs 
also decrease with increasing height, specifically in 77 percent of the 
cases up to 700 mb. Such vertical gradients of baroclinicity in the low 
troposphere agree with the generally-held concept of frontal zones over 



continental areas. 
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Table 10. Statistics for the numerically-analyzed HBZs in 
Figs. 32-34. All numbers represent averages. 







< 


— Level (mb)-- 


- > 


Frontal Parameter 




1000 


850 


700 


500 


width of HBZ, An (deg lat) 




6.2 


6.6 


6.8 


7.6 


A0* through HBZ (C) 




22 


16 


13 


12 


A0*/An (C [deg lat] ) 




3.5 


2.4 


1.9 


1.6 






< - 


— Layer (mb)-- 


> 






1000/850 


850/700 




700/500 


A0* along HBZ boundary (C) 


K: 


11.5 


9.6 




11.6 




W: 


5.0 


6.8 




10.0 


A0*/ AZ along HBZ boundary 












(C [1000 ft]-l) 


K: 


2.3 


1.9 




1.5 




W: 


1.0 


1.4 




1.3 






< 




> 






1000/850 


850/700 


1000/700 


Slope of HBZ boundary 


K: 


1/55 


1/175 




1/115 




W: 


1/15 


1/164 




1/90 



Thermal Character of HBZ Boundaries --HBZ boundaries are known to 
show an associated increase of potential temperature with increasing 
height. O’Connor (United States Navy, 1952) observed rather small 
changes, as 298 to 302 K along the wintertime polar front from sea level 
to 500 mb, respectively. However, the Canadians indicate larger changes 
are likely along frontal zones (Godson, 1950), as 296 to 312 K in winter 
and 305 to 323 K in summer along the polar front from 850 to 500 mb, 
respectively. The 0* values used here steadily and consistently increase 
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in magnitude along the HBZ boundaries at a rate which averages close 
to 1.3 C (1000 ft) from 850 to 500 mb, and hence are quite compatible 
with the Canadian values. The less conservative nature of the cold 
boundaries is borne out by the higher average rate of compared to 

that of the warm boundaries in the layer 1000 to 700 mb (i.e., 2.1 versus 
1.2 C [1000 

Frontal Slope --The frontal-zone slopes are greater at low levels 
and for the warm-side boundary of HBZs , as expected from conventional 
notions of fronts. However, the slopes appear more erratic than those 
which would be expected from manually -drawn fronts; negative slopes (i.e., 
HBZs sloping toward warm air with increasing height) occur in 20 percent 
of the cases. In the layer below the 700 mb level, such behavior is 
partly due to the poor numerical resolution inherent in the mass structure 
model for the layer 775 to 1000 mb (Chapter III) and the predominance of 
steeply-sloping cold fronts in the sample considered here. In the layer 
above the 700 mb level, the picture may be further complicated by fronts 
which are rooted in the high troposphere and have position and intensity 
not intimately associated with the low-level phenomena (Fig. 35). A 
similar feature was evident in the DG0* analyses for the summer-season 
case study. 

Time Continuity --The cross-sections of HBZ-boundary profiles 
shown in Figs. 32-34 are graphed in such a way as to aid in analyzing 
time continuity. Table 11 presents information on two time-continuity 
parameters derived from changes of the warm (w)- and cold (k)-side 
boundaries of the HBZs in the layers 1000/850 mb and 850/700 mb. 

1. The percentage of time that the algebraic sign of the 12-hour 
slope change is the same for the w and k boundaries of the HBZ. 
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Fig. 35. Schematic profiles of HBZ boundaries for coexisting 
low- and high-troposphere frontal zones (Endlich 6e McLean, 1957). 
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Table 11. Time-continuity statistics for numerically-analyzed 
HBZs in Figs. 32-34. Column a--Algebraic signs of 12-hour slope 
change same for w- and k-side boundaries of HBZ . Column b--Algebraic 
signs of 12-hour change in HBZ width and slope of positive and negative 
DG0* axes differ. 



Cross section group 


Average 


frequency (percent) 




a 


b 








+ axis 


-axis 


I 


100 


67 


67 


II 


50 


67 


100 


III 


50 


33 


33 


IV 


90 


64 


53 


V 


100 


70 


33 


VI 


100 


50 


50 



2. The percentage of time that the algebraic signs of the 12-hour 
change in HBZ width and slope of each of the DG0* axes differ. 

The internal consistency of the numerical product is shown by the 
fact that over 80 percent of the time slope changes in the positive and 
associated negative axes of the DG0* field are similar in algebraic sign. 
There is no apparent pattern in the changes of slope and width of HBZs. 

In other words, as the frontal-zone boundaries steepened or became more 
horizontal with time the HBZ widths may either increase or decrease. 

The relation of 12 -hour changes in HBZ width and magnitude of 
DG0^ is quite conclusive in the case of increasing HBZ widths, since 73 
and 87 percent of the DG0^ values along the associated positive and nega- 
tive axes, respectively, decreased during the same period. However, the 
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situation for decreasing HBZ width is more nebulous with an average of 



Thus, the weakening HBZs (using zone width as an indicator) correlate 
well with changes in the magnitude of DG0*, but the cases of intensifying 
HBZs are indefinite in this respect. This situation is not unexpected as 
a consequence of poor resolution in the numerically-processed data 
(Chapter II). 



continuity (0000 GMT 26 January to 1200 GMT 27 January 1967) of 850 mb 
numerically-analyzed frontal positions. Only those frontal phenomena 
over the ample-data United States area are in focus here. The aim is 
one of illustration, not statistical comparison. The time continuity of 
the HBZ warm-side boundaries appears logical in the area east of the Rocky 
Mountains. However, behavior over the mountainous southwest United States 



only 46 percent of the changes in 




being positive for both axes. 



The final figure of this chapter (Fig. 36) depicts a two-day time 



is erratic. Further, the absence of occluded structure is to be noted 
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Fig. 36. Time continuity of 850 mb numerically-analyzed fronts 
(i.e., warm-side boundaries of HBZs): 0000 GMT 26 January 1967 (thin 

solid lines); 1200 GMT 26 January 1967 (thin dashed lines); 0000 GMT 
27 January 1967 (heavy solid lines); 1200 GMT 27 January 1967 (heavy 
dashed lines) . 



CHAPTER VI 



THE DG0* MODEL: 

A CASE STUDY OF HYPERBAROCLINIC ZONES 
IN RELATION TO CLOUD AND VERTICAL MOTION 

The hyperbaroclinic zones of major interest to the field 
meteorologist are those migratory transition regions of the middle 
latitudes separating cold polar from warm tropical air. Such zones 
are invariably associated with synoptic-scale cloud and vertical motion 
patterns which have a common signature in the large-scale sense but vary 
in small-scale details mostly due to geographic influences. The clouds 
are the visible manifestation of saturation due to relative moisture 
content and the time-integrated effect of vertical motion. The 
intensity and shape of the vertical motion field is largely a function 
of the two baroclinicity-de pendent parameters which may be described as 
the vertical rate of change of vorticity advection and the horizontal 
Laplacian of thickness advection (Krishnamurti , et al., 1966; Sanders, 
1963). However, a one to one correspondence between baroclinicity and 
cloud patterns or baroclinicity and vertical motion patterns does not 
exist due to the diverse combinations of the thermal and vorticity forcing 
functions. Nevertheless, the relation of cloud to hyperbaroclinicity 
has been modeled for use in synoptic analysis and forecasting, starting 
with the Norwegian concept of polar front- phenomena in the 1920s 
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(Petterssen 1956a, b) and more recently by those observant meteorologists 
attempting to exploit the operational uses of the weather satellite 
(Boucher and Newcomb, 1962; Anderson et al., 1966). 

Fig. 37 qualitatively implies a cloud-HBZ relation widely held 
by the current-day field meteorologist. Though certainly less ideal 
than the ear ly-century concepts, the model must be described as a 
heavily-smoothed version of the real atmosphere. The extensive cloud 
masses and bands made up of cirriform, stratiform, and cumuliform clouds 
near the cloud -mass vortex, in the warm sector, and especially adjacent 
to and within the frontal zone, are indicators of significant tropos- 
pheric vertical ascent. For other areas, especially the closed-cellular 
and cloud-free regions, the vertical motion is likely to be negative 
over at least part of the troposphere. 

Until the existence of a model to quantitatively specify the 
location and intensity of frontal zones, an objective comparison between 
hyperbaroc linicity , vertical motion, and cloud has not been possible. 

The research described below concerned itself with such a substudy, with 
emphasis on the relation of satellite-observed cloud to the HBZ . 

Selection of Period and Area of Study 

An eleven-day winter period, 7-17 December 1966, was chosen to 
investigate the relationships of satellite-observed clouds and FNWC’s 
vertical motion fields^ (Hughes, 1967) to the HBZs derived from the 

^FNWC calculates vertical motion from a diagnostic (jo equation in 
the form v2(ocD)+(pfr|/R)(SWSp2) =V^(g/f J[D,T])-J(T,n) - J(D,g/f [V^T]), 
where a is a stability parameter, J is the Jacobian, D is the isobaric- 
height deviation from the United States Standard Atmosphere, r| is absolute 
vorticity, and other symbols are standard. 
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DG0* model using the 0* appropriate to the layer 1000/700 mb. The 
weather satellite information used in the study consisted of NMC- 
analyzed nephana lyses of the AVCS-type ESSA III observations. In 
order to consider a significantly large area while at the same time 
maintaining reasonable concurrency of cloud and HBZ analyses, only 
the nephana lyses within four hours of the 0000 GMT and 1200 GMT times 
were utilized. 

Use of the satellite-observation periods of 2000 to 0400 GMT 
and 0800 to 1600 GMT determined two major geographical areas of study. 

The time interval centered at 0000 GMT gave coverage over most of the 
North Pacific Ocean plus a portion of North America west of 95W, while 
the second interval gave coverage of the eastern North Atlantic Ocean 
east of 55W , Europe, and western Asia to 90E. Since the seasonally- 
imposed northern limit to the daylight satellite view is close to 60N 
and since the numerical frontal analysis is not too reliable at latitudes 
less than 15N due to sparse data, the south and north boundaries of the 
areas of study were chosen accordingly. 

The December 1966 weather situation was characterized by 
abnormally strong zonal flow over much of the Northern Hemisphere 
with mean troughs near both coasts of North America (Dickson, 1967). 

In the North Pacific area a series of frontal-wave cyclones moved 
through the latitude band 40-55N, recurving northeast into the Gulf of 
Alaska. In the North Atlantic region, the major storm track was con- 
fined to polar latitudes north of 60N due to the presence of a relatively 
strong, stationary subtropical high in the middle latitudes. Because of 
the blocking situation, only the trailing portions of fronts crossed 
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into the North Atlantic area of study. Over the Eurasian continent, 
the weather in the early part of the period was dominated by a high 
pressure cell over Siberia and weak frontal lows in the Mediterranean. 

In the last half of the period the weather was generated mainly by low 
cells in the North Sea and over the northern part of the continent with 
frontal activity across Europe extending southward to the Mediterranean. 

In the illustrations and cloud statistics which follow only 
cloud amounts are considered. That is, no attempt was made, on the 
basis of the nephanalyses , to make any distinction between high, middle, 
and low cloud types. An exception to this procedure occurred for the 
cases of cloud bands or masses composed entirely of cirriform cloud. 

Such clouds were considered to have no direct relation to the low- 
troposphere HBZs and were not included in the statistics. With reference 
to standard satellite meteorology terminology, the cloud amount covered, 
C, indicates > 80 percent cloud cover while mostly covered, MCO, specifies 
50-80 percent sky coverage. 

Patterns of Vertical Motion and Cloud in Relation 
to the Numerically -Analyzed Fronts 

Some notion of the empirical relations between the numerically- 
analyzed 1000 mb front, satellite-observed clouds, and 850 mb vertical 
motion may be obtained by viewing a 3-day sequence (15-17 December 1966) 
for the North Pacific area (Figs. 38-40). The area is not necessarily 
the most desirable for such a demonstration due to sparsity of data, but 
the wide expanse of uniform surface aids in generating patterns of a 



more universal nature. 
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Cloud banding appears to be well-related to the 1000 mb fronts 
in orientation. However, as expected, the colocation is not as ideal 
as that of the model in Fig. 37, The low- latitude fronts are associ- 
ated with long-trailing cloud bands extending from middle latitudes 
into the subtropics. Many of the HBZ boundaries are within the area 
covered by cloud masses which include bands not evident to the cloud 
analyst. It is to be noted that not every DG0* axis is simply related 

to a cloud mass or band nor does every such cloud entity correspond to 

- -2 

a synoptic-scale HBZ, The vast majority of DG0* maxima > 15 x 10 

(100 km) are found in the covered and mostly covered areas. 

The relation of vertical motion to HBZs is more difficult to 

analyze. As discussed earlier, a given pattern of vertical motion is 

related more clearly to cloud than to the zones of hyperbaroclinicity , 

2 

However, major areas of 850 mb vertical motion (i,e. centers > 1,0 
cm sec are somewhat related to HBZs, with 60 percent of the centers, 
of ascent located within such zones and most of the remainder in the 
adjacent warm-air subbaroclinic zone, similar to models suggested by 
Sanders (1963) and Todson (1964), The orientations of low- latitude verti- 
cal motion patterns do not appear to be closely associated with either the 
cloud bands or DG0* axes, A part of the problem of interpretation here is 
due to lack of radiosonde data in the area. 

In summary, the illustrations shown here do not suggest an 
unqualified set of usable relations between HBZs, vertical motions. 



2 “1 
This area was contained within the 1 cm sec isoline for 500 mb 

vertical motion more than 95 percent of the time. Patterns at the two 

levels were very similar; only the magnitude differed. 
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and clouds. However, it is obvious that HBZs and areas of 850 mb 
vertical motion equal to and greater than +1.0 cm sec ^ are associ- 
ated with major cloud bands and masses. Further, there is the 
suggestion that the less significant (i.e., very narrow) frontal cloud 
bands are found adjacent to and oriented similarly to the positive DG0* 
axes; in these regions the vertical motion is generally less than + 1 
cm sec In order to quantitatively specify the relations between 

cloud and HBZ , the next section presents statistics derived from the 
7-17 December 1966 period. 

Satellite-Observed Cloud Bands Versus Hyperbaroclinic Zones 
Stratification by Latitude 

Statistics on cloud bands, stratified by latitude, are shown 
in Table 12. Average DG0* values are indicated for each of the three 
latitude bands in the two major hemispheric regions sampled, namely, 
the North Atlantic east to western Asia and the North Pacific east to 
western North America. The relatively high DG0* intensities along the 
positive axes indicate that the warm-side boundaries of hyperbaroclinic 
zones represent closer approximations to a first-order discontinuity 
than the cold-side boundaries, as concluded in Chapter II. Also, the 
average magnitude of the DG0* values increases with latitude, an in- 
direct measure of the poleward intensification of the frontal zones. 

The statistics on clouds versus HBZs demonstrate that the 
average percentage of C and MCO cloudiness in the HBZ is larger than 
that in the SBZ at corresponding latitudes, consideri*ng all latitudes 
between 15N and 60N. Conditions in the low latitudes of the Pacific 
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area and the high latitudes of the Atlantic area are most responsible 
for the near two to one ratio of cloud percentages in the HBZ versus 
the SBZ . Although the combined percentage of C and MCO cloud cover 
inside (especially) and outside the HBZ markedly increases with lati- 
tude, the actual area covered by clouds in each latitude band increases 
more slowly since the area ratios of the two latitude bands, 30-35N 
and 45-60N, to 15-30N are .86 and .65, respectively. 



Table 12--Relation of sate llite -observed clouds to hyperbaro- 
clinic and s ub baroc linic zones, stratified by latitude, for period 
7-17 December 1966. 



Latitude Average DG0* Percentage of 
band (deg) (C [ 100 km]“2) latitude band 

within HBZ 



Percentage 
of HBZ with 
cloud amount 



Percentage 
of SBZ with 
cloud amount 





Negative 

Axes 


Positive 

Axes 




C 


MCO 


C 


MCO 






Atlantic 


-Euras ian- 


-African Region 






15-30N 


.11 


-.08 


39 


13 


12 


8 


12 


30-45N 


.19 


-.09 


30 


41 


18 


19 


15 


45-60N 


.20 


-.19 


37 


65 


9 


30 


7 




Area -we 
average 
15-60N 


ighted 

for 


37 


31 


13 


16 


12 






Pacific 


-North American Region 








15-30N 


.11 


-.05 


24 


30 


16 


11 


8 


30-45N 


.17 


-.13 


37 


51 


26 


25 


16 


45-60N 


.20 


-.15 


33 


63 


15 


38 


14 



20 19 



Area-weighted 
average for 
15-60N 



30 



41 



11 
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Stratification by 
Geographical Area 

Table 13 considers each type of cloud band (frontal covered, 
non-frontal covered, and non-frontal mostly covered) in each of the 
following six geographical zones: A, North Pacific; B, western North 

America; C, eastern North Atlantic; D, Europe; E, western Asia (C, D, E 
poleward of 30N) ; and F, At lantic-Africa-Asia , 15-30N, Here, the word 
^’frontal” refers to the satellite cloud entity labeled as such by the 
ESSA nephanalyst, apparently when the band is categorized as C and is 
thought to be associated with a conventional front as depicted on the 
NMC surface chart. 

Pacific area A indicates that the average angle between the C 
bands and their associated DG0* axes is larger (29 deg versus 15 deg) 
and has a wider range (10-90 deg versus 5-45 deg) in the frontal than 
in the non-frontal cases. As discussed in Chapter III, the operational 
numerical frontal analysis model does not adequately depict occlusions 
while the amplitude and orientation of frontal cloud bands relates well 
to such fronts; this is a possible explanation for the observed 
relation. Corresponding values in area C appear to disagree with this, 
since the average angle in the frontal-covered category is small. 
However, as noted earlier in this chapter, the storm tracks and their 
attendant occlusions were well within the area studied in the Pacific 
but not in the Atlantic. Therefore the Atlantic and Pacific statistics 
do not conflict. 

Further examination of Table 13 discloses that in the frontal- 
covered category a greater portion of the cloud band fell into the HBZ 
than in the SBZ in both ocean areas A and C, but the cloud amount in the 
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Table 13--Relation of satellite-observed cloud bands to hyper- 
baroclinic and s u bbaroclinic zones in the latitude band 15-60N, 
stratified by geographical area: North Pacific--A; North America, 

west of 95W--B; North Atlantic, poleward of 30N, east of 55W--C; 
Europe, poleward of 30N--D; Asia, west of 90E--E; and the Atlantic- 
Africa- Asia area, 15-30N between 55W and 90E--F, 



Area Type cloud Number 
band and of 

coverage bands 

sampled 


Average 
angle 
between + 
DG0* axis 
and cloud 
band (deg) 


Range 

of 

angles 

(deg) 


Average DG0* 
(C[100 km]-2) 

along 

positive axis 


Percentage 
of cloud 
bands in 
HBZ 




Frontal 


C 


18 


29 


10 to 90 


.17 


55 


A 


Non-frontal 


C 


9 


15 


5 to 45 


.13 


35 




Non-frontal 


MCO 


8 


8 


0 to 10 


.09 


33 




Frontal 


C 


1 


90 




.20 


40 


B 


Non-frontal 


c 


0 












Non-frontal 


MCO 


0 












Frontal 


c 


8 


14 


5 to 35 


.25 


61 


C 


Non-frontal 


c 


10 


29 


0 to 90 


.16 


51 




Non-frontal 


MCO 


1 


10 




.07 


50 




Frontal 


C 


2 


12 


10 to 15 


.12 


78 


D 


Non-frontal 


C 


7 


13 


0 to 30 


.16 


63 




Non-frontal 


MCO 


0 












Frontal 


c 


0 










E 


Non-frontal 


c 


10 


28 


0 to 90 


.23 


55 




Non-frontal 


MCO 


2 


12 


10 to 15 


.15 


20 




Frontal 


C 


0 










F 


Non-frontal 


C 


0 












Non-frontal 


MCO 


12 


Indefinite 




.11 


50 
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SBZ is by no means insignificant. This result agrees with the inter- 
pretation of the model in Fig. 37, the broadening effect on the 
numerically-analyzed fronts, and the extensive zone width due to using 
0* rather than 0*. 

In the non-frontal C and MCO classes combined the statistics 
are inconclusive. In areas A and E most of the cloud banding occurs 
in the subbaroclinic region while in D this is true of the hyperbar o- 
clinic zone; in C and F the bands are evenly split between the HBZs and 
SBZs . Non-frontal cloud bands, especially MCO, are frequently associ- 
ated with the stratus and stratocumulus attending the stable regions of 
anticyclones and, as such, large amounts in the SBZs are to be expected. 
The satellite cloud bands associated with major positive DG0* 

axes were also studied. An axis was classified as major if at least 
— -2 -2 

one DG0* value > 15 x 10 C (100 km) occurred along it. The results 
were similar to those presented in Table 13 and are not shown here. 



CHAPTER VII 



THE DG0* MODEL: HEMISPHERIC FREQUENCY AND AVERAGE INTENSITY 

OF HYPERBAROCLINIC ZONES 

The DG0* analyses yield quite readily to an objective calcu- 
lation of the frequency and average intensity of HBZs . Two 30-day 
periods were selected to demonstrate a synthesis of objective 1000 mb 
frontal analyses for the Northern Hemisphere; one, 15 January through 
14 February 1966, exhibits winter-time conditions, and the other, the 
month of July, 1967, shows mean summer-time conditions. For both winter 
and summer periods the frequencies of subjectively-analyzed NMC surface 
fronts were also calculated. This was done over a reduced area com- 
prising the western Northern Hemisphere, as appropriate to NMC*s 
facsimilie transmissions of the 0000 GMT and 1200 GMT surface and 
1000/500 mb thickness fields. 

The frequency and intensity statistics were accumulated manually 
from the 0000 GMT and 1200 GMT numerical DG0* analyses (1000/700 mb model) 

by tabulating the occurrence of and average intensity along those positive 

— -2 

axes where DG0* > .05 C (100 km) for each of 3844 grid squares. Each 
square measures 381 km on a side, true at 60N. The NMC frontal positions, 
not allowing an intensity calculation, were tabulated in a similar way for 
frequency only. 
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Winter-season Sample: 15 January through 14 February 1966 

Fig. 41 highlights the important features of frequency and in- 
tensity of numerically-analyzed 1000 mb fronts for the period selected. 

For frequency, the maximum possible occurrence per grid square is 62 
which is equivalent to the number of 0000 GMT and 1200 GMT synoptic times 
involved in the statistic. Minimal frequency is zero, indicating no 

occurrence of numerical fronts in the grid zone during the period. Aver- 

-2 

age intensity values have a lower limit of .05 C (100 km) since only 

i 

those positive DG0* axes with values equal to or greater than five units 
were counted. The upper limit of intensity is only a practical one since 
most map-time intensities are less than 60 units of DG0*. In general, 
analysis at latitudes less than 20 deg lat is omitted due to the unreli- 
able nature of data for these sparse-data areas. 

As background to the interpretation which follows, it is to be 
noted that the last half of January, 1966, was anomalously cold over North 
America while the first half of February, 1966 (that is, at least from 
February 5 forward), represented a trend back to normal winter-time con- 
ditions (Ludlum, 1966). This suggests that over the ample-data North 
American region frontal activity should show supernormal occurrence at 
relatively low latitudes with subnormal, but not zero, frequency at polar 
latitudes. In fact, this same trend was evident over all of the western 
North Hemisphere (Posey, 1966; Green, 1966). The eastern section of the 
hemisphere, especially Europe and western Asia, was associated with rela- 
tively more disorganization and meridional flow as compared to the western 
section, particularly in January. The numerically-derived frontal averages 
for these months show essentially the same features. 
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Numerical Frequency 

Axes of maximum frequency are multiple with mostly zonal 
orientation, except along the continent-ocean interface where overtones 
of coastal orientation are evident. Frequencies equal to or in excess 
of 30 percent are concentrated over land areas. Regional maxima appear 
on the Canadian west coast near 55N, over the extreme North Pacific in 
the Kamchatka Peninsula area, east of Japan at 35N, over the western and 
southern Mexican coastal area, and along the Asiatic front in the foothills 
of the Himalayas. 

There are several areas of zero or near-zero frequency, as noted 
from the 3-percent lines in Fig. 41. The more extensive areas of minimal 
frequency, as in various parts of the Arctic latitudes, may be regarded 
as the traditional air-mass source regions. However, certain localized 
areas of this type are more difficult to explain, especially in view of 
the high frequencies in the surrounding grid zones. Examples may be found 
over the east coast of the United States near 35N, southwest United States, 
southern Finland, and ocean areas in the latitudinal band 45-50N. Most of 
the local minima over land have the common feature of being in a mountain- 
lee area while such ocean areas more likely represent tenuous boundaries 
between arctic air excursions toward the south and maritime- polar air 
traverses toward the north. Minimal zones south of the most equatorward 
positive DG0* axes are not indicated. These zones are generally equator- 
ward of 20N and display erratic patterns due to sparse data. 

Numerical Intensity 

The intensity field shows that the intense maritime fronts occur 
farther poleward than their counterparts over the continental regions. 
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Fig. 41. Frequency and intensity of numerically-analyzed 1000 
mb fronts in the period 15 January through 14 February 1966. Frequency 
isolines for 3^ (2 cases), 20^ fl2 cases), 30^ (19 cases), and 45^ (28 
cases). Frequency axes: 20-29^, heavy dashed lines; > 30^, heavy 

s oJid lines. Areas with frequencies < 3^, shaded. Axes of average 
DG0* intensity (units, 10“^ C [100 km]“^): 15-20 units, light dashed 

lines; > 20 units, light solid lines with crosses. Intensity areas > 
30 units enclosed by dashed line. Intensity axes omitted through areas 
of frequency < 3^. 
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In general, highest intensities are found over land with absolute 

-2 

maxima of .49 and .47 C (100 km) on the Alaskan west coast near 60N 
and 55N, respectively. Such figures are of significance considering 
that a total of 11 and 23 cases, respectively, are associated with these 
intensities . 

The large-scale relation between intensity and frequency is one 
of maximum intensity lagging the most frequent position of the front. 

That is, the former is toward the coldest air from the axis of greatest 
frequency. 

NMC Frequency 

It is of interest to compare the frequencies of NMC’s manually- 
analyzed fronts (Fig. 42) to the objective version in Fig. 41. The format 
in both figures is the same for easy comparison. However, an intensity 
comparison is not possible in the case of the NMC fronts. 

The major differences between NMC and numerical in this short- 
period mean are similar to those found in the single-time views. For 
instance, the many occlusions in the eastern ocean areas give more 
meridional orientation to the NMC compared to the numerical frequency 
axes since the latter does not portray such fronts very well (Chapter 
III). Another extensive area of difference is eastern Canada and the 
North Atlantic Ocean, particularly Greenland and Hudson Bay, which shows 
high frequency in the DG9* field, but near zero occurrence on the NMC 
chart. Such a difference is to be expected since it is not common practice 
to depict orographic and other semi -permanent frontal zones on NMC's daily 
operational weather charts. The NMC and numerical patterns are most com- 
patible over the Canadian Archipelago, mountainous sections of western 
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Canada, and the United States and the North Atlantic equatorward of 
40N. 



Summer -seas on Sample: July, 1967 

Data for the summer-season month, July 1967 (Fig. 43) , are pre- 
sented in the same fashion as the winter period discussed above. 
Initially, the data were stratified by the upper-air observation times 
0000 GMT and 1200 GMT in order to analyze diurnal variations, but one 
month’s data proved to be inadequate for this purpose. 

To aid in interpreting the July 1967 frontal statistics, the 
following comments regarding the Northern Hemisphere weather are of 
interest. The subject month was much like January 1966 in that North 
American temperatures were well below average east of the Rockies. 

Record lows were common and associated with polar-front penetrations 
well equatorward of normal, especially over southeast United States 
(Ludlum, 1967). In the European and central and west Atlantic and 
Pacific areas, low-troposphere high pressure systems were well developed 
with ant icyc Ionic and supernormal zonal flow at higher latitudes than 
usual for July. Thus, frontal positions may be expected poleward of 
normal at these longitudes. In the Asiatic sectors the flow patterns 
were not suggestive of any marked or unusual frontal phenomena (Dickson, 
1967 b). 

Numerical Frequency 

Compared to winter, there are fewer and less extensive frequency 
axes with values > 20 percent but relatively more points with frequency 
> 30 percent along the axes. Isolines of 45 percent and 60 percent are 
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Fig. 43. Frequency and intensity of numerically-analyzed 1000 
mb fronts, July 1967. Freauency isolines for 3^ (2 cases), 20^ (12 
cases), 30*^ (19 cases), 45^ (28 cases), and 60^ (37 cases). Frequency 
axes: 20-29^, heavy dashed lines ; > 305 S, heavy solid lines. Areas with 

frequencies < 3^ shaded. Axes of average DG0* intensity (units, 10*2 
C [100 km]“2). 12.5-15 units, light dashed lines; > 15 units, light 

solid lines with crosses. Intensity areas > 20 units enclosed by 
dashed line. Intensity axes omitted through areas of frequency < 3^. 
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to be found over Asia, Africa, Mediterranean Europe, and mountainous 
North America. An absolute maximum frequency of 70 percent is reached 
along the southern periphery of the Iranian mountains (near 30N) , just 
inland from the Persian Gulf and the hot Saudi- Arabian desert to the 
southwest . 

Areas with near-zero frequencies are of greater extent in summer 
than in winter, both in the Arctic and as buffer zones between the 
zonally-oriented major frontal axes. The multiplicity of such minimal 
areas points out the reduced latitudinal migration of fronts in the 
former as compared to the latter season. 

Except for northern sections, frontal frequencies of any con- 
sequence over oceans are non-existent. However, the relation of 
latitudinal positions of continental and maritime fronts is similar 
to winter. The circulation patterns prevailing in July 1967 undoubtedly 
emphasize this relation. The zonal extent of fronts is interrupted by 
the coastal orientation and orographic influence over the western contin- 
ental areas of Europe and North America. The extensive Asiatic front 
has shifted to the polar-side desert basins of the Himalayas at most 
longitudes, a seasonal shift expected from the nature of the thermal 
patterns in the two seasons . 

The frontal ring over mountainous western North America may be 
explained by a number of features on the periphery of a July maximum- 
temperature area normally centered near 37N, 113W (U. S. Navy, 1966). 

One such feature was discussed in Chapter V (Fig. 26), namely the dif- 
ferential heating across the land-sea interface. Other sections of the 
ring- like front are due to the polar front, usually found near the 
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Canadian-United States boundary in the summer season, the sparse-data 
and cooler mountain and ocean sectors along the south-southwest fringe 
of the desert area, and a combination of moisture and temperature differ- 
ences between land and water areas to the east and southeast of the hot 
elevated basin. The las t -enumerated phenomenon generates the dew-point 
front, occurring often through the Texas area in summer (U.S. Navy, 1952), 
A less prominent but similar situation occurs over northwest Africa. 

Numerical Intensity 

As expected, the intensity field is more complex and much reduced 

in magnitude compared to winter. The maximum average intensity of .21 C 
-2 

(100 km) , found just off the Atlantic coast of African Morocco, is less 
than half of the absolute maximum noted in the winter sample. However, 
intensities over mountainous western North America closely approach 20 
units. Most zonally-oriented fronts at latitudes less than 40N and most 
ocean fronts are not associated with intensity axes of significance, using 
> 15 units as the criterion. The lag of the maximum intensity zones 
relative to maximum frequency zones, evident in winter over the east con- 
tinent-west ocean areas, is not apparent in the summer picture; rather, 
the orientation of frequency and intensity axes tends toward coincidence 
in the mean. 

N>KI Frequency 

The pattern of axes on the NMC analyses are not entirely according 
to expectation (Fig. 44). The multiplicity of major axes of frontal 
frequency not only exceed NMC*s winter picture but numerical in summer as 
well. The relation to numerical is least expected as a result of vari- 
ations noted in the daily patterns (Chapter III), An area in point is 
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eastern North America and extreme western Atlantic, usually one of 
maximum agreement in the NMC and numerical frontal analyses. Since 
the manual product does not suffer greatly from sparse-data boundary 
effects and since the seasonal stationary fronts are not depicted in 
the daily subjective analyses, the frontal frequency pattern over south- 
western North America is simpler than numerical. The usually trouble- 
some extreme northeast Pacific area yields different solutions on the 
subjective and objective versions. In general, central and western ocean 
areas and the Arctic region show the best comparison between NMC and 
numerical patterns , 

Prior to the development of the frontal analysis model described 
here, a simple, objective, quantitative method for measuring the location 
and intensity of fronts did not exist. Thus, frontal climatology has 
remained qualitative and descriptive over the years. However, it is 
obvious from the results of the prototype experiment of this chapter 
that quantitative frontal climatology is now feasible and that such a 
synthesis may be expected to relate well to geography, season, and the 
daily weather sequences. 



CHAPTER VIII 



RESUME* AND PROPOSALS FOR FURTHER STUDY 



Summary of Results 

The following major accomplishments were achieved during the re- 
search study reported on here. 

1. The most important result of the investigation is the develop- 
ment of the first operationally successful numerical-objective model for 
the synoptic-scale analysis of the location and intensity of hyperbaro- 
clinic zones (HBZs). The scheme locates the constant pressure HBZ 
boundaries (i.e., fronts) along the axes of maximum baroc linicity gradi- 
ent, simulating, in principle, the positioning of fronts by subjective 
analysts who adhere to the transition-zone concept of frontal phenomena. 

The model is easily adapted to a computer-oriented national weather ser- 
vice operation, as evidenced by the 1000 mb Northern Hemisphere numerical 
frontal analyses and prognoses being produced twice daily by the Fleet 
Numerical Weather Central, Monterey, California (FNWC) . 

2. A four-dimensional synoptic and statistical diagnosis of HBZ 
structure and motion is made possible by the development in 1 . For 
this purpose, modeled and real-atmosphere virtual potential temperature 
data were employed to illustrate the relation of the merits and deficiencies 
of . the . analysis model to the numerically-processed and computed data. Gen- 
eral characteristics of the HBZs are shown to be universal, although a 
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particular mode of numerical-data processing and computation, namely 
that of FNWC, was used throughout the study. Grid size , frontal- 
zone intensity (especially that due to zone width), air-mass 
dimension (i.e., width of the adjacent subbaroclinic zone), and use of 
layer versus mandatory- leve 1 temperatures are shown to be factors in 
yielding computed frontal structure which differs in degree from the 
real atmosphere as a function of the numerical manipulation of data. In 
general, the HBZs are weakened and the widths broadened by numerical pro- 
cessing, with these effects somewhat proportional to the original intensity 
in the transition zone. For a given mesh size, use of mean-layer vice 
single- level temperatures reduces the broadening effect on the HBZs. In 
this regard, use of the 1000/700 mb mean virtual potential temperature - 
field, 0*, was found to be suitable for generating 1000 mb (or surface) 
fronts. The reduction in mesh length, which must be closely tied to data 
density, appears to offer the most potential for improvement in accuracy 
and detail of numerical frontal analysis. 

3. Several case studies from the years 1964-67 in extended time 
and space frames indicate the broader utility of the HBZ as compared to x 
the conventional air-mass associated front. The numerical- ana lysis model 
has the advantage of locating fronts at any level, whether moving or oro- 
graphically fixed, and not just those fronts separating air masses of 
recognizable thermal characteristics and dimensions. However, occluded 
structure is not evident from use of the model without empirical modi- : 
fication. Further, the relation of HBZs to vertical motion and cloud is 
shown to be complex and not operationally usable for locating HBZs. 

4. The study also proposes a more elegant model for the purpose 

of analyzing the mesoscale structure of HBZs. The high resolution inherent 
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in vertical soundings, combined with a vertical extrapolation procedure, 
increases the resolution and amount of baroclinicity data to the extent 
that the fine structure of baroclinicity is clearly evident on the con- 
stant pressure mapping. The capacity of the major national weather 
centrals is not sufficient presently to take full advantage of such a model 
but it is offered as a prototype for future operational use. 



Proposed Areas for 
Further Research 

The following areas for extended research on HBZs appears pro- 
ductive at this time. 

1. It has been demonstrated that use of vertically-resolved 
thermal data allows the mesoscale structure of HBZs to be analyzed* In 
reality, the synoptic-scale HBZs comprise a multiple number of these sub- 
zones of hyperbaroclinicity (Kreitzberg, 1964)* In this study such res- 
olution was introduced at the baroclinicity stage of the HBZ analysis. 

It is proposed that the more basic field, temperature, be resolved by 
extrapolating the vertical temperature structure into the horizontal or 
constant pressure plane and then applying the subject model for the 
numerical analysis of HBZs* With a significant increase in data density, 
a judicious grid-size reduction can be fully effective in increasing the 
resolution of the analyzed HBZs* For a high-resolution model of this 
type, it is presumed that considerable smoothing of temperature and/or 
its derivatives will be necessary in order to locate the boundaries of the 
synoptic-scale HBZs, as found by previous researchers working with detailed 
frontal structure (e*g*, Beniura, 1962)* 
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2. The operational surface-front analysis model proposed here, 
utilizing 0*, is essentially a thermal- as opposed to a weather -oriented 
scheme* With little guidance from national or international meteorologi- 
cal commissions, the latter approach is becoming popular with the field 
meteorologist as a convenient way of forcing apparent agreement with the 
satellite-observed cloud masses and bands. Perhaps regarding the frontal 
zone as the centrum of a circulation entity, (Bleeker, 1961; Eliassen, 

1962; Todson, 1964; and Kirk, 1966) would allow relating vertical motion, 
and hence cloud, more intimately to baroclinicity , thus joining the thermal 
and weather concepts of the frontal zone. With the advent of night and day 
coverage by a multiple number of weather satellites, empirical modeling of 
vertical circulation, weather, and baroclinicity may now be feasible. 

3. The limited substudy concerned with the long-period (as 30 
days) frequencies and average intensities of HBZs indicated the potential 
use of these charts for discriminating between the fronts fixed to oro- 
graphy and those which are progressive. For operational analysis and fore- 
casting it may be desirable to eliminate the orographic-type fronts in 

the daily analysis of HBZs. Some preliminary regional work contributing 
to the solution of this problem has been done by Reed (1961), who deter- 
mined tables of average seasonal surface temperatures for synoptic map 
times over the western United States. Reed suggests, but does not show, 
that the deviation of temperature from a long-term mean be used to locate 
fronts in the complex mountain areas. The existence of monthly mean 
radiosonde temperatures, recently updated by the United States Navy (1966), 
makes such an experiment feasible. It is proposed that the HBZ model 
developed here be used in connection with the temperature anomalies. 
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4, Again, with reference to the substudy mentioned in 3 above, 
a computation of the average frequency and intensity of numerically- 
analyzed HBZs , as determined from a suitable period of several years, 
will allow the establishment of a much needed quantitative climatology of 
Northern Hemisphere fronts. 

5. Operationally, there is need for a computer-plotted front in 
addition to the isolines of the front- location parameter since there is 
some subjectivity in locating the implied boundaries of the HBZ on 
meteorological charts of small scale. Some work has been done on this 
problem, using a directional third derivative of 0*, but the effort was 
terminated due to the errors inherent in fields of such high-order 
derivatives . 

Final Remarks 

Experience in applying the numerically -analyzed fronts emphasizes 
the need for a man-machine mix (Sne liman, 1969) due to the current gross 
grid size and poor horizontal resolution in the numerically-processed 
operational data. Using the DG0* model, the computed frontal positions 
may be regarded as reasonably good first approximations to the actual 
locations, Mesoscale adjustment of the objective ly- located fronts as a 
function of the adjacent and updated pressure, wind, and temperature 
observations is still deemed to be necessary. Such an approach does 
require plotted data as well as a knowledgable meteorologist but the 
time and effort involved is much less than a totally manual approach. 

It is of interest to note the recent theoretical work by Williams 
(1967) which has as its goal the dynamical modeling of mesoscale fronto- 
genesis from initial large-scale conditions of a simple, vertically-bounded. 
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atmosphere. Williams* developmental model yields vertically- 
sloping frontal zones whose boundaries may be located in terms of 
the characteristics of potential temperature and wind profiles, quite 
analogous to the operational frontal-analysis model described here. 



APPENDIX A 



PRELIMINARY EXPERIMENTS IN HY PERBAROC LIN IC -ZONE ANALYSIS 

Experiments Using Wet-bulb Potential Temperature 

The initial impetus for the study of numerical hyperbaroclinic 
zones (HBZs) came from some excellent statistical work on the air-mass 
distributions of wet-bulb potential temperature, 0^, published by 
Harley (1962). 0^ has been the prime parameter in the daily opera- 
tional frontal-contour analyses emanating from the Central Analysis 
Office, Meteorological Service of Canada (CAO) , since 1953. 

Harley presents, among other things, monthly means, standard 
deviations, and modal values of 0^ in three air masses described by the 
Canadians as maritime tropical, maritime polar, and maritime arctic. His 
data sample is provided by a careful analysis of a full year (1957) of 
radiosonde observations from 89 stations in the area 25-90N and 15W to 
170E across North America. The method used to determine an air-mass 0 

w 

allows his statistics to be interpreted as usable climatological indicators 
of the warm-air edge of the HBZs separating two adjacent but dissimilar air 
masses. As an example. Table 14 contains the monthly means and standard 
deviations of frontal 0^s for the polar front, which separates maritime 
tropical from maritime polar air. The 0^ values are taken from Harley, 
but rounded to the nearest .1C. 

Some preliminary work, attempting to assess the utility of 0^ as 
a front- location parameter, is described in the following subsection. In 
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order to utilize the empirical results of Harley and yet limit the initial 
experiments on locating HBZ boundaries, the following conditions and pro- 
cedures were imposed: (1) 0 and V0 are considered as the prime front- 

w w 

location parameters, (2) only the poleward boundary of tropical air (i,e,, 
polar front) is located, (3) the 850 mb level is a suitable single test 
level, and (4) since objectivity is the goal, any successful subjective 
analysis scheme should be a prototype for a numerical-objective approach. 

Table 14. Monthly mean 0 and 0 and their standard deviations, 
for the polar front, as analyzed Sy CAO.^ 



Month 


®w 

(°C) 


Standard 

deviation 

(°C) 


0 

e 

(°c) 


Standard 

deviation 

(°C) 


Jan 


14.0 


1.1 


42.1 


3.1 


Feb 


14.0 


1.0 


42.2 


3.0 


Mar 


14.2 


1.0 


42.6 


2.9 


Apr 


14.9 


1.0 


45.0 


3.0 


May 


15.5 


1.3 


46.7 


4.0 


Jun 


16.8 


1.3 


50.8 


4.4 


Jul 


17.7 


1.0 


53.9 


3.7 


Aug 


17.7 


1.1 


53.9 


3.8 


Sep 


17.4 


0.9 


52.9 


3.0 


Oct 


16.5 


1.0 


49.8 


3.3 


Nov 


15.6 


1.0 


46.8 


3.3 


Dec 


14.9 


1.1 


44.8 


3.3 
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Relations of 0 and V0 to the Polar Front 
w w 

Summer Case Study: 

Manual Analysis 

The National Meteorological Center (NMC) and CAO 850 mb polar- 
front analyses over North America, for 35 synoptic times (0000 GMT 21-26 
May 1963, 0000 GMT 21 June to 1200 GMT 26 June 1963, and 0000 GMT 12 July 
to 1200 GMT 16 July 1963), were screened for their relationship to the 
0 field. The 0 value at each radiosonde station was obtained graphi- 
cally from the United States Naval Weather Service’s standard thermo- 
dynamic chart, the AROWAGRAM, using transmitted radiosonde temperature 
and dew-point data. 

Specifically, answers to the following equations were being 
sought from analyses of the 850 mb polar front on the NMC and CAO charts: 

1. What is the relation of NMC and CAO 0 s along the front? 

w 

2. Are values of 0^ relatively uniform along either of the 

fronts ? 

3. Are frontal values of 0 conservative with time? 

w 

4. In view of answers to _1, _2, and do Harley’s data assist 
in objectively placing an air-mass boundary? 

With reference to these questions, analysis of the May, June, and 
July 1963 data indicated the following conclusions. The numbered para- 
graphs refer to the foregoing questions. 

1. A sampling of the polar-front 0^s at longitude intervals of 
5 deg over North America generated the percentage frequencies and averages 
shown in Table 15. Similarity of NMC and CAO frontal 0^s does exist in 
the mean, but apparently there are differences in 0^ values along the two 
versions of the polar front on a given day. 
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Table 15. Percentage frequency of wet-bulb potential temperature, 
0^, along the 850 mb polar fronts analyzed by NMC and CAO. 



0^ Percentage Frequency 



(°C) 


21-26 


May 


1963 


21-26 


June 


1963 


12-16 July 


1963 




NMC 




CAO 


NMC 




CAO 


NMC 


CAO 


12 








1 




3 






13 


















14 


7 






5 




2 






15 


3 






11 




7 






16 


13 




7 


14 




10 


4 




17 


15 




28 


12 




12 


16 


22 


18 


36 




36 


22 




46 


30 


28 


19 


23 




26 


14 




9 


26 


26 


20 


3 






16 




9 


13 


7 


21 






3 


2 




2 


4 


2 


22 








3 






7 


11 


23 
















4 


Average 


17.5 




17.5 


17.4 




17.4 


18.6 


18. ( 


Number of 
Cases 




39 






57 




46 





2. Inspection of the average frontal values of 0^ as a function 
of longitude, by analysis center, shows the geographical variation of 0^ 
along the front (Table 16). The variation is greatest in the case of NMC 
fronts. Apparently the polar-front values of 0^ are not uniform on the 
average, as indicated in Table 15 as well. 



Table 16. Average wet-bulb potential temperature, 0 , as a function of longitude, along 
the 850 mb polar fronts, analyzed by NMC and GAO. 
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3, Some measure of the conservativeness of frontal 0 values 

w 

and other features may be obtained by noting a characteristic day in 
the series. Fig, 45 for 0000 GMT and 1200 GMT 21 June 1963 shows the 
relation of NMC and CAO polar fronts to 0^ data and Harley’s statistical 
estimate of the frontal positions. 

4, Fig, 45 shows that the orientation of manually-analyzed and 
statistical frontal positions are similar but not coincident. In general, 
the temporal non-conservativeness and spatial non-uniformity of 0^ sug- 
gests the probable failure of trying to objectively locate fronts based 

on a single climatological value of 0^, But the preponderance of NMC and 

Canadian positions within one standard deviation (a) of the monthly mean 

0 is suggestive of a statistical approach in which the warm-air boundary 
w 

of the HBZ is searched for within the area bounded by 0 + k, where 

k < 3 C, 

Winter Case Study: 

Numerical Analysis 

The 1963 summer tests were considered at least encouraging enough 
to continue experimenting with the application of statistical values of 0^ 
for frontal analysis. However, since the aim was to develop an objective 
system operationally applicable to the Northern Hemisphere, use of numeri- 
cally-processed data and fields^ became necessary and computations of 0^ 
were converted to the numerical mode. In this connection an analytical 
expression for 0^ (Berry et al, , 1945) of the following form was programmed 
for solution on FNWC’s Control Data Corporation 1604 digital computer. 



Data were processed in accordance with FNWC's mass structure model as 
described in Chapter III. 
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Fig. 45. 850 mb polar front analysis. 1200 GMT 21 June 1963: NMC front, heavy solid line, 
CAO front, thin solid line. Modal position (Harley, 1962), dashed line; 9+1 standard deviation 
(Table 14), dotted lines. NMC front, 0000 GMT 21 June 1963: VAV. 1200 GMT NMC front at same 0^ 
as the 0000 GMT position; VAT, 
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0 =0 exp - i - 

w e V 



^-[ihsLIq^ Vs 1000' 



0 C 

w p 



1000 C 



D. 



where 



. X m 

P c 



(25) 



(26) 



m is mixing ratio, L is latent heat of condensation, subscripts s and c 
indicate saturation and condensation levels, respectively, e is vapor 
pressure, 1000 is 1000 mb, and other notations are standard. Due to the 
limited nature of this experiment the various complex numerical approxi- 
mations and iteration schemes used to compute 0^ will not be discussed 
here. Such are on file at FNWC. 

The synoptic time 1200 GMT 5 February 1964 was chosen for experi- 
mentation in objective frontal analysis using the numerically- produced 
wet-bulb potential temperature field. In this case the field is actually 
wet-bulb virtual potential temperature, 0*, since FNWC's processed tempera- 
tures approximate virtual temperatures. Fig. 46 shows a portion of the 
850 mb 0* field over North America. The area within 0^ + 3 q, appropriate 
to February (Table 14), is regarded as the most likely region within which 
the polar front is located. Using 0* as a suitable approximation to 0^, 
the first approximation to the position of the polar front is the zone 

(called search zone) between the 0* isolines of 11, OC and 17. OC, Follow- 

w 

ing the definition of the front used here (Chapter I), the line of maximum 

V0* within the search zone became the second approximation to the polar 

front. In this case, the orientation of the axis of maximum V0* is simi- 

^ w 



lar to the subjectively-located fronts, especially CAO, but locations differ. 
This and other tests of a like nature, although suggestive of the direction 
in which to pursue the research, did bring about the abandonment of 
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wet-bulb potential temperature as a front locator for the following 
reasons , 

1. 0* was being considered solely because of the existent 

w 

air-mass statistics of Harley (1962). By using the Smithsonian 

Meteorological Tables (List, 1963) the 0^ statistics may be readily 

converted to their 0^ counterparts. Table 14 includes the monthly 

mean frontal values, 0 , so determined. 

e 



due to the multiple numbers of scans necessary to approximate 0* to an 
acceptable degree of reliability. Consequently, the next testing phase, 
still following the lead of the Canadian statistical work, utilized 
numerically-computed 0* vice 0*. 

Experiments Using Equivalent Potential Temperature 

For reasons discussed above, further tests with 0* were aborted 

w 

in favor of the more easily computed equivalent virtual potential tempera- 
ture, 0*, using a form of the analytical expression in Eq. (26). 



2, 0* computations are difficult and costly in computer time 



One of the front- locator parameters developed is of the form: 




i 



where i refers to £olar, maritime, or arctic fronts; lv0*| represents 
the magnitude of the gradient of 0*; and 0*^ is a normalized value of 0*, 
defined as 



9 *" = 0 * 



e* - 9* . 

e e . 



(28) 



e. e 



i 



e . 

1 



For the purpose of this experiment, 0* + ka were regarded as equivalent 

_ ^i 

to 0 + ka. The location of the front is taken as the axis of maximal 



values of F. . 

1 
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For the purpose of applying the scheme of Eq, (27) as a third 

approximation to the 850 mb polar front of 1200 GMT 5 February 1964 

(Fig. 46), the following procedure is implemented. First, the 0* lines 

e 

are redesignated as 0*^ lines. Thus 0* of 33. 2C, 42. 2C, and 51. 2C 

e e 

become 0^^ of 33. 2C, 42. 2C, and 33. 2C, respectively. Next, a search 
e 

zone defined by 0^ + 3a is scanned for the axis of maximal values of 

e — 

F^ . Such a procedure forces location of the front between the maximum 

V0* and 0*, however always within the search zone. To assure this con- 
e e 

dition in the machine computation, the F^ values may be arbitrarily set 

to small values outside the search zone. In those cases in which 0* does 

e 

not lie in a direction toward warm air from the maximum V0* the F. com- 

e 1 

putation of the front location is less accurate than that dependent on 

the maximum V0* alone. 

e 

Enhancement factors --In order to ensure placement of the front to 
the warm-air side of the axis of maximum V0*1 in the search zone, a 
weighting factor, W, of dual form was used in combination with F^ to gen- 
erate the final front-location parameter considered in this statistical 

synoptic approach. Thus, 



F' . = F. W = 0*" |ve*| w, 






(29) 



and, as for F., the axis of maximal values of F’. locates the polar front. 
’ 1 1 ^ 



At the boundary of and inside the search zone defined by 0* + 3a 



W = 




(30) 



e 

where X is a constant for a particular front. Thus, increases from 
one at the cold-side boundary to (0^+3a)/ (0*“3a) at the warm-side boundary 
of the search zone. X values up to nine were used to enhance the search- 

area gradients of 0^ 

e . 
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Outside the search zone. 



W = W2 




(31) 



e 

which always yields a fraction. X values up to twenty were used to reduce 

the gradients outside the search area to extremely small fractional values 

As an example, following this scheme in the 1200 GMT 5 February 

1964 case, Table 17 shows the weight factors and other information which 

result for the specified 0* lines. For cases where maximum V0* is at the 

location 0* + 3a, F*. does not change the front- location estimate. For 
e 1 

all other locations of a V0* maximum, F*. yields a more consistent front- 

e ’ 1 

position estimate than F^ or maximum V0* alone. Fig. 46 shows the fourth 
approximation to the polar front using the F\ approach. 



Table 17. Table of related values from Eqs . (27) and (29), as 
applied to a February 850 mb polar front. 





0* 

e 

(°C) 


0**^ 

e 

(°C) 


Symbolic weight 
factor 


Magnitude of 
weight factor 


F* , 
1 




24.2 


24.2 


^2 


20 


0.31 


0 -3 
e 


33.2 


33.2 




9 


1.00 


Q 

e 


42.2 


42.2 




9 


1.31 


0 +3 
e 


51.2 


33.2 




9 


1.69 




60.2 


24.2 


^2 


20 


0.31 



Empirical adjustment to the statistical approach is almost endless 
Wind shear or other front-associated parameters may be introduced as com- 
ponents of the weighting factor in Eq . (29), various values of x niay be 
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employed in the enhancement scheme, and so on. A definite problem is 
proper hemispheric combination after each of the three fronts, polar, 
maritime, and arctic, is determined separately. Various arithmetic 
schemes of combining the individual or fields were tried. The 

results were not satisfactory. 

Since a more elegant, as well as a more utilitarian, model was 
developed (Chapter II), somewhat forced by the experiences just recounted, 
the illustrations here will be truncated to one experimental numerical 
plot of the arctic frontal field (Fig. 47), using a form of the enhance- 
ment scheme just outlined (Eq . 29). The axes of the elongated F’ isolines 

a 

represent the arctic front. Note that this approach yields hemis pherically 

continuous and mostly zonally-oriented fronts. Comparison with the 

Canadian arctic front is good, as might be expected from the dependence 

of F* on Harley's values of 0*; these values are used to some extent by 
a e 

the Canadians in their frontal analysis approach. A more extensive 
analysis of Fig. 47 is not warranted in view of the resume’ immediately 
following . 



Summary 

The various frontal-analysis schemes described in this Appendix 
restrict the fronts to a definite number, as well as to specific zones 
determined by air-mass statistics (Harley, 1962). That is, fronts out- 
side the zone of 0^ + 3a polar, maritime, and arctic fronts are not 
allowed. This is a very artifical restriction. Thus, final placement of 
a front is not dependent solely on factors involved in the definition of 
a front accepted here (see Chapter I) but rather on air-mass labels. A 
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Fig, 47. Frontal field, F^ , generated for use in locating the 
850 mb arctic front at 0000 GMT 5 March 1964. The search zone boundaries 
given by F^ = 1: heavy lines. The axis of F^ within search zone locates 

the estimated position of the arctic front. The CAO version of the arctic 
front is shown for comparison: xxx . 
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dynamic-synoptic approach to the problem is desirable. Such is the 
model finally adopted, as described in Chapter II, 

Further, 0* and 0* were abandoned as front- location parameters 
w e 

in favor of 0* due to the relative inaccuracy of numerical moisture 
analyses. It was found that motorboating values and/or very low dew 
points were improperly rejected in the numerical processing procedures. 
This has the effect of giving excessive 0* and 0* values in the relatively 
cold dry air, thereby reducing the 0* and 0* gradients and causing the 
misplacement of the frontal positions. Additionally, detail inherent in 
the moisture field was not resolved by analysis procedures tuned to the 
gross mesh size appropriate to current-day numerical analyses. It might 
be added that some hope of improved and detailed moisture analysis has 
been given recently through at least partial incorporation of satellite 
data in the numerically-processed fields. However, objective frontal 
analysis experiments with parameters dependent on detailed moisture 
information are considered premature even to this date. 
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